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Significant progress has been achieved under grant AFOSR-86-0031 in understand¬ 
ing the physics underlying ultrafast transient phenomena that occurs in the sem¬ 
iconductor microstructures. This understanding is essential for making the neces¬ 
sary advances for the future generation of ultrafast microelectronic devices. 

Microstructure samples have been obtained primarily from H. Morkoc of University 
of Illinois and other groups such as M. Niigaki of Hamamatsu Photonic KK,TRidl32s 
Bertaska of McDonnel Douglas, and Emil Koteles of GTE. 

The following summarizes our accomplishments in the major effort on III-V sem¬ 
iconductor microstructures and a minor effort on Il-VI compound. 

1. Nonequilibrium Phonon Effects on the Energy Relaxation and Lifetime 
of Photogenerated Carriers in GaAs MQW 

We have carried out time-resolved photoluminescence with the finest time- 
resolution of 2ps performed to date to investigate the ultrafast processes which 
occurs in undoped GaAs/AlGaAs microstructures. 

The existence of a large population of nonequilibrium (hot) phonons after an initial 
rapid carrier cooling in an undoped multiple GaAs quantum well (55A) structure 
excited by 500fs laser pulses is experimentally verified by time-resolved photo- 
luminescence studies using a 2ps-time-resolution streak camera. The hot phonon 
population is directly monitored by measuring time-resolved phonon replica 
luminescence below n=I electron to heavy-hole transition energy. It reaches to a 
maximum in ~ 3ps after the end of laser pulse and decays with a 30 ps time con¬ 
stant. Its effects on carrier dynamics are the following two aspects. First, the 
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energy relaxation of carriers after an initial rapid cooling is substantially 
suppressed due to re-absorption of hot phonons by carriers. Second, the carrier 
density decay is strongly influenced by the hot phonons in the system which 
stimulate phonon replica emission. This stimulated phonon replica emission is 
mainly radiated and highly polarized in the well plane (x-y plane). An effective 
carrier depletion time [density decreases by a factor of c -1 from n e (t =0)] is 
obtained to be as short as lOps. 


This work is of particular importance to the semiconductor community because it 
solves the long standing controversial issue regarding the phonon assisled transi¬ 
tions below n=l electron to hole transition in GaAs quantum wells that was ini¬ 
tially suggested by Professor Holonyak’s group (see for example, Phys. Today, 40, 
Oct. 1980). We have explicitly demonstrated that the hot phonons emitted by hot 
electrons can stimulate phonon-replica-emission. 


Part of this work was accepted to be published in Phys. Rev. B. 


2. Observation of Ultrafast Lateral Diffusion of Photogenerated Carriers in 
An Asymmetric GaAs Quantum Well 


An experimental investigation of in-plane expansion of phologenerated carriers in a 
single asymmetric GaAs quantum well is carried out using a streak camera and a 
subpicosecond laser source at 4.3K. It is found that the carrier expansion process is 
enhanced in the well plane by restricting carrier expansion in the well direction. 
The effective diffusivity D' of photoexcited carrier in the asymmetric well has 
been directly determined to be 10 ( 'an 2 l s at 4.3K, which is about four orders 
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magnitude larger than the thermodynamic value in bulk GaAs. The expansion 
velocity is determined to be 3.2xl0 8 cm I s , which is four times larger than the 
initial Fermi-velocity of the electrons. 

The observation of an anomalously fast carrier expansion process in an asymmetric 
quantum well, but not in the symmetric well is most important. It provides a 
new insight to the problem of how to design semiconductor structures in 
order to obtain the desired carrier transport. 
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This work was submitted to Phvs. Rev. Ixtt. for publication. 

3. Determination of Band Offsets via Optical transitions in Ullrathin 
Quantum Wells 

This work demonstrates what experimental data from optical transition should be 
used and how it is sensitive to the Q -value of band offset. We have shown that 
the precise Q -value can be determined by systematically measuring the n=l light- 
and heavy-hole energy spacing (A£ ) in the sensitive zone of well width (15 to 
80A). Furthermore, it is also shown that different connection rules and hole 
masses yield different Q -value and the appropriate connection rule can be discrim¬ 
inated also by systematically measuring A£ as well as the energy spacing between 
ground and first electron eigen state (AE 12 ) in the sensitive zone of L.. 

This work was published in Phys. Rev. B, 37, 1048 (1988). 

Recently, we are investigating the influence of different connection rules on the 
exciton binding energy in III-V semiconductor quantum wells, and hence on the 
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determination of band offset based on the data of A E. 

4. Physics in semiconductor GaAs and GaSe under laser-driven shock-wave 
compression 

Using a pump-probe technique, the photoluminescence peak of bulk GaAs was 
observed to blue shift and split into two components under the laser-induced shock 
wave loading condition. The two components correspond to the transitions from 
the T 6 conduction band to the valence heavy- and light-hole subbands, because of 
symmetry breaking by the uniaxial shock compression along the [001] direction. 
From the blue shift, we deduced our picosecond-laser-induced shock pressure of ~ 
10k bar. 

In GaSe, we observed a broadening of the spontaneous emission which is propor¬ 
tional to the shock pressure and attributed to a shock-wave-induced exciton colli¬ 
sion mechanism due to the directional motion of particles in the shocked region. 

These works were published in Appl. Phys. Lett. 51, 1789 (1987) and Appl. Phys. 
Lett 52, 93 (1987). 

5. Optical transition and recombination lifetime in quasi-zero dimensional 
electron system in CdS x Sc 

We report on the observation of optical transitions between quantized levels 
(IS ,\P ) in the conduction and valence bands in quasi-zero dimensional electron 
systems in CdS x Sc by 4K and 300K photoluminescence measurements. From 
the relative luminescence intensities arising from the IS — IS and IP —1 P 


■r>, 










§ 

s 

i 

$ 

i 

£ 

i 


transitions, it is shown that the photoexcited electrons (holes) occupied in the 
different quantized levels in conduction (valence) band are not coupled, or only 
weakly coupled. Picosecond luminescence studies of the optical transitions using a 
streak camera reveal a significant influence of three-dimensional confinement on the 
recombination probabilities of photogenerated electrons and holes. The recombina¬ 
tion life time of electrons and holes in quasi-zero dimensional system decrease with 
the diameter (d) of system from 210ps for d=10.2nm to 70ps for d=7.4nm at 4K. 
The ratio of recombination life time for IS —IS transition and IP —IP transition 
seems to be independent of the diameter of system, and is measured to be about 
3.5. This ratio can be interpreted by the calculated squared transition matrix ele¬ 
ments and the nonradiative processes. 

This work was published in Appl. Phys. Lett. 51, 1839 (1987). 
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Effects of nonequilibrium phonons on the energy 
relaxation and recombination lifetime of photogenerated 
carriers in undoped GaAs quantum wells 
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Kai Shum, M. R. Junnarkar, H. S. Chao, and R. R. Alfano 

Institute for L'ltrafast Spectroscopy and Lasers 
Departments of Electrical Engineering and Physics 
The City College of New York, \ew York, New York, 10031 

and 

H. Morkoc 

Coordinated Science Laboratory, University of Illinois at 
Urbana-Champaign, Urbana, Illinois 61801 


ABSTRACT 

Time-resolved and time-integrated photoluminescence studies of 
an undoped multiple GaAs quantum well structure excited by 0.5ps 
laser pulses have revealed several important experimental observations 
on the behavior of photoexcited carriers. A large population of none¬ 
quilibrium longitudinal optical phonons produced in the energy relaxa¬ 
tion process of hot carriers manifests itself by the nonequilibrium pho¬ 
non stimulated phonon-replica which located at ~30meV below n=l 
electron-hole transition. The energy relaxation is substantially 
suppressed due to the existence of nonequilibrium phonons after an ini¬ 
tial rapid cooling (0-5ps). The number of photoexcited carriers 
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decreases anomalously fast within the first JOps after the excitation by 
laser pulse. An effective carrier depletion time is determined to be as 
short as lOps. A mechanism which leads to such a short carrier deple¬ 
tion time is associated with the nonequilibrium phonon stimulated 
phonon-replica. 

PACS numbers: 63.20.Kr, 72.20.Jv, 73.40.Lq. 
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I Introduction 

Recently, there have been growing interest in studying the dynamics of photo¬ 
generated carriers in semiconductor quantum well structures since it reflects the fun¬ 
damental interactions among electrons, holes and phonons. These interactions deter¬ 
mine the performance of ultrahigh speed electronic and optoelectronic devices. In the 
case of bulk semiconductors, the slowing of hot carrier relaxation is attributed to 
either screening 1 of electron-phonon interaction or re-heating of carriers by nonequili¬ 
brium (NE) phonons produced in the relaxation process of hot carriers.”’ 3 A compara¬ 
tive study 4 of the hot electron cooling rates in undoped multiple quantum well 
(MQW) structures and in bulk GaAs by time-resolved measurements of optical 
absorption and gain concluded that the rates were approximately same at a photogen¬ 
erated carrier density of 2.5xl0 17 cm' 3 . This is expected from a simple theory 5 in 
which lattice is treated as a heat bath for quasi-equilibrium carriers. In the case of 
modulation doped MQW structures, quasi-steady-state experiments have been carried 
out from two groups 6,7 to investigate the interaction of electrons and holes with pho¬ 
nons. These experiments have generated conflicting results regarding the presence 6 or 
absence 7 of NE phonons. Time-resolved PL measurements with ~20ps time resolution 
were reported by Ryan et al . 8 who found that the cooling of hot carriers was 
anomalously slow after 40ps. 

Reviewing the aforementioned previous works two important points were 
noticed. First, one can not extract information on the initial carrier relaxation process 
either from quasi-steady-state experiments or from time-resolved measurements with 
time resolution greater than lOps. To assess the importance of NE phonon effect on the 
energy relaxation of hot carriers one uses a corresponding theory to match the experi¬ 
mental carrier temperature cooling curve or the logarithm of power loss as function 
of the inverse of the carrier temperature. There has been no report on the direct evi¬ 
dence for existence of NE phonons in quantum well structures under high photoexci- 
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tation. It is necessary to study the initial carrier relaxation process to substantiate the 
existence of NE phonons and its importance on the hot-carrier relaxation since the 
phonon lifetime is about 7ps. Second, one can not easily obtain information on photo- 
generated carrier lifetime in modulation doped MQW structures because the large 
doped-in carrier density is comparable to the photoexcited carrier density. Further¬ 
more, carrier lifetime can influence the cooling of hot carriers. 9 Therefore, it is neces¬ 
sary to study photogenerated carrier dynamics in undoped MQW structures in order 
to obtain information both on the energy loss rates and on the carrier lifetime. 

In this paper, experimental observations are reported from the measurements of 
time-resolved photoluminescence (PL) with a 2ps time resolution and time-integrated 
luminescence spectra from an undoped GaAs MQW structure with dependences on 
lattice temperature, excitation intensity, and polarization. The NE phonons emitted 
by hot electrons are directly observed by measuring the time-integrated as well as 
time-resolved phonon replica luminescence lying below the n=l electron to heavy- 
hole transition energy. The energy relaxation of hot electrons is found to be substan¬ 
tially suppressed when a large population of NE phonons is present after an initial 
rapid cooling. The photoexcited carrier density extracted from a fitting of time- 
resolved PL profiles at different emitted photon energies decreases non-exponentially 
and very rapidly within the first 30ps. An effective carrier depletion time is deter¬ 
mined to be as short as lOps. The mechanism leading to such short carrier depletion 
time is associated with NE phonon stimulated phonon-replica. 

11 Sample 
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The undoped GaAs/Alj.Gaj.^.As MQW structure investigated was grown by 
molecular-beam epitaxy on a (001 )-oriented undoped GaAs substrate. The MQW con¬ 
sists of 50 periods of 55A thick GaAs and 100A thick Al 03 Ga Q7 As layers, and fol¬ 
lowed by a 1.2 um GaAs buffer layer. The lateral size of the structure is about 
2x4mm. The good quality of sample is confirmed by well resolved heavy-hole and 
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light-hole excitonic structures of room temperature PL spectra. The FWHM for the 
heavy-hole excitons is about 13meV. The sample was mounted on the cold finger in 
an optical helium cryostat. 

Ill Experiment 

An ultrashort laser pulse of 0.5ps duration at 620nm was used to excite the 
electron-hole pairs in the sample. Most of light is absorbed in both the GaAs wells 
and the AlGaAs barriers. Photogenerated carriers in the barriers will either diffuse to 
the GaAs buffer layer or be captured into the wells. The latter will result in uncer¬ 
tainty about carrier density in the wells by a fact of ~3. Residual transmitted light 
is absorbed in the GaAs buffer layer. The diameter of excitation area is about 160wm. 
To avoid damaging of the sample surface, the maximum excitation power is adjusted 
such that photogenerated carrier density is on order of 10 19 cm' 3 . The 0.5ps laser 
pulse was generated from a colliding-pulse passive mode-locked dye laser and 
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amplified by a four-stage dye amplifier pumped with a frequency doubled Nd: YAG 
laser operated at 20Hz. In order to keep the time resolution of Hamamatsu streak 
camera system within 2ps, PL was spectrally resolved using different narrow band 
filters. The luminescence intensities were corrected for the nonlinearity of streak 
rate, the spectral response, and the transmission of each narrow band. For the time- 
integrated PL measurements, a Stanford boxcar and a GaAs photomultiplier were 


IV Results and Discussion 
1. Steady-state 

Time-integrated luminescence spectra from the MQW structure excited by the 
0.5ps laser pulse at various lattice temperatures (7 L ) are shown in Fig.l. These spec¬ 
tra were taken in a conventional backward Raman configuration r (y ,y )z , where z is 
the growth direction. In order to eliminate radiation from the edges of the sample. 
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the luminescence spot was first imaged on an aperture with aid of a streak camera 
and then re-focused onto a vertical slit of.a grating spectrometer. Several features are 
displayed in the data shown in Fig.l. (1) The emission band A on high energy side of 
the spectra arises from the recombination of photogenerated n=l electron-hole. The 
peak of A shifts towards the low energy side as the T L increases. An high energy 
tail on the peak A develops with increasing T L . 

(2) A broad emission band C at low energy side of the spectra arises from the 
GaAs buffer layer. The total emission intensity of C decreases as T L increases. 

(3) The most interesting feature of the spectral data is the appearance of an 
emission band B, located about 30meV below the n=l electron-heavy-hole transition. 
This B emission band is attributed to the NE phonon stimulated phonon-replica. The 
following four reasons support this assignment to the B band: 

(i) This B emission band does not appear in the PL spectra taken at low power 
excitation about lW/cm 2 at 4.3K using a cw 488nm line of an argon-ion laser. How¬ 
ever, the PL spectra using a very weak train of laser pulses (l20fs) directly from the 
colliding-pulse passive mode-locked dye laser as the excitation source, with an excita¬ 
tion power density in range of 10~ 5 —1CT 3 W / cm 2 , shows a weak electron to acceptor 
emission band separated by 17meV from the n=l electron-heavy-hole transition 
accompanied by its phonon-replica 10 at low temperatures. This extrinsic emission 
band disappeared completely when the sample temperature was raised to ~ 80K, 
whereas the B band displayed in Fig.l exists up to room temperature. Moreover, the 
relative time-integrated PL intensity of the B to A band under 0.5ps light pulse exci¬ 
tation shown in the Fig.l decreases as the excitation power density decreases. The B 
band in Fig.l is hardly visible when the excitation power density is lower than 
10~ 3 /’ m , where P„, —10 12 W / cm 2 is the maximum value of the excitation power den¬ 
sity. This contradicts what is generally expected for impurity emission. Since the 
concentration of acceptors is low in our sample as confirmed by the luminescence stu- 
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dies using the weak train of 120fs light pulses, an impurity emission band should be 
more readily apparent and pronounced at lower excitation. Therefore, the B emission 
band shown in Fig. 1 can not be attributed to electron to acceptor luminescence under 
high excitation. 

(ii) The intensity ratio between the B and A emission bands as indicated by tri¬ 
angle in Fig.2 increases as lattice temperature increases. To convince the emission 
band B arises from NE phonon stimulated recombination process, we calculate the 
sum of occupation numbers both for the equilibrium phonons ( lattice temperature 
Ti ) and for the NE phonons ( carrier temperature T c ) by assuming the carrier tem¬ 
perature is same as the effective temperature for the NE phonons after 30ps. This 
will be discussed later with the carrier temperature. The calculated result is shown 
by the solid line in Fig.2. The fitting of the total phonon occupation number to the 
obtained intensity ratio from the spectra is impressively good. This implies that the 
intensity of band B is well correlated with the LO-phonon population. Because of the 
participation of NE phonon the intensity of B band (see also time-resolved lumines¬ 
cence at 780nm in Fig.4) is strongly enhanced especially at low temperatures. This 
further supports our assignment of the B band to the NE phonon stimulated phonon 
replica. 

(ili) It is also not possible to attribute the emission band B to the emission from 
renormalized band-band transition because the peak position of the A band (766.7nm) 
and the B band at 4.3K do not change with variation of excitation intensity by a fac¬ 
tor of ~ 10 3 . A further support to above statement is the fact that the spectral posi¬ 
tion of peak A exactly coincides with each other for the two different luminescence 
studies using different light excitation sources: one is the weak 120fs-pulse train with 
a repetition rate of 125MHz, and the other is the amplified 0.5ps-pulse with a 20Mz 
repetition rate. 

(iv) We have also studied the polarization of the emission bands A, B. and C by 
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measuring "right-angle" 11,1 “ z (y , x )y and z (y , z )y integrated luminescence spectra. 
It was found that the B band emitted at the sample edge along the y-direction (see 
inset of Fig.3) was highly polarized in the x-direction. The A and C emission bands 
were depolarized, independent of excitation power density and lattice temperature. 
The intensity ratio between the x- and z-direction for the B band is about 20 in the 


1 

W 

& 

1 


lattice temperature range from 4.3 to 300K at full excitation power density P ,„. But, 
the ratio is strongly dependent on the excitation power density at a given lattice tem¬ 
perature. For example, at T L = 20A' , the ratio decreases from 20 to —6 as the excita¬ 
tion power density decreases from P,„ to 0.017P,,,. Fig.3 shows the spectra for 
T l = 100A' at two different excitation intensities. The lower solid trace was taken 
in z (y , x )y configuration (TE polarization) at the excitation density of 0.0l7/’ m . It is 
identical to the spectrum taken in z (y, z )y configuration ( I'M polarization) which is 
not shown. The B band disappears due to the absence of large population of NH pho¬ 
nons at the low excitation and high lattice temperature. The upper solid and broken 
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traces are the TM and TE spectra at the full excitation P n ,, respectively. The inten¬ 
sity ratio of TE and TM remains 20. excitation P m . It should be pointed out that the 
energy position of A band in Fig.3 is lower than that in Fig.l by ~4meV while the 
position of B band remains unchanged. This low energy shift of A band is due to the 
self-absorption effect. The polarization behavior of B band is consistent with the 
results from Raman scattering experiments reported by Zucker et al. ] ~ where the 
pump photons are provided by the external laser source. In our case the pump pho¬ 


tons arise from the recombination of n=l electron and n=l heavy-holc at the subband 


edges. This process is strongly enhanced by the presence of a large number of NE 
phonons emitted by hot-electrons at high excitation. The depolarized feature of the 
band C is expected since the emission is from the GaAs buffer layer. However, the 
polarization behavior of the emission band A is not understood. In terms of the selec¬ 
tion rule of dipole recombination described by Iwamnra ct u/. 13 , the emission intensity 








of n=l electron-light-hole for the TE polarization should be four times larger than 
that for the TM polarization. There should be no TM emission arising from n=l 
electron-heavy-hole recombination. It should be mentioned here that the B band in 
the Fig.l which was detected along the z-direction is depolarized and much weaker, 
at least a fact of 20, than detected along the y-direction. The former is expected due 
to the symmetry and the latter provides an important rule for understanding a fast 
carrier density decreasing process which will be discussed in the timc-rcsolved section. 

In order to further substantiate our assignment of the B band, two additional 
questions must be addressed: 

Cl) Why no phonon replica appears on the high energy side of the A band, if a 
large NE phonon population really arises? A high energy replica would arise from a 
recombination of an energetically elevated electron (hole) by absorbing a phonon and 
hole (electron). However, the elevated electron (hole) will be very quickly scattered 
by the other hot electrons (holes) through a strong electron-electron, electron-hole, 
and hole-hole interaction before it recombines with a hole (an electron). Therefore, 
the high energy replica can not be observed in PL spectra. But, this re-heating process 
should result in a retardation of hot carrier cooling which is indeed consistent with 
the time-resolved PL data described in the next section. 

(2) Why is the energy separation between the peaks of A and B bands (/; Ag ), on 
order of ~ 30meV which is smaller than expected (30meV) 7 The 30m eV energy 
separation in fact gives a further strong support of our assignment of the B band. In 
an equilibrium state, the interaction strength of electron and LO phonon in a bulk 

GaAs is proportional to where Q is the phonon wave vector. Therefore, the pho¬ 
nons in vicinity of Q = 0 would be expected to more strongly couple with electrons 
giving rise to a phonon replica at just 36meV below its primary emission band. In 
our highly photoexcited MQW sample, LO phonons are driven to a NE state by a rapid 
initial energy relaxation process of photoexcited hot electrons. The develop of these 
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NE both in time domain and in two-dimensional wave vector space (</*, parallel to the 
well plane) has been recently studied by several groups. 14,15 Two important results 
on the number of NE LO phonon population are : (i) At given q , it can reach to a 
maximum in about 0.5 - 3ps and then decrease, (ii) At a given time, it increases stee¬ 
ply from a minimum of q min , reaches to a maximum at q c , and then decreases. What 
is essential concerning the question is that there is a maximum NE LO pho¬ 
non population which is located both in a narrow wave vector space of the 
vicinity of q c and in the well plane and lasts in a short time period. This 
portion of NE phonons behaves like coherent Bosons. The existence of the "coherent 
Bosons" will tend to increase the number of the Bosons (near q c ) in the system with 
a rate proportional to the present number of Bosons resulting in a stimulated phonon 
emission process. This process can give rise to a phonon replica emissions, when an 
electron in the conduction subband recombines with a hole in the valence subband, 
not only a photon but also a phonon or more phonons will be emitted to join the 
Boson system. In order to conserve the momentum of electron-phonon system in this 
recombination, the electron must possess a momentum of q r (neglecting the niomen- 

ft ~q e ~ 

turn of the emitted photon) which is associated with a kinetic energy of — 

Therefore, the energy E AB should be equal to LO phonon energy minus E which is 
about 6meV corresponding to the values of 1.03xl0 6 cm _l and 0.067m „ for q r and elec¬ 
tron effective mass m,. Thus, this gives 30meV for E AB . 

Over years, there has been much debate 13,16 " 18 over the interpretation of the 
spectral features locate below the n-1 electron to heavy-hole transition energy, 
mainly because many species exist in this energy region: LO-phonon replica, impurity 
states, and many body band gap renormalization. Holonyak and co-workers 17 have 
demonstrated the LO-phonon participation in GaAs-AIGaAs based QW laser emission 
by observing mre than one LO-phonon sidebands in the laser operation below the 
n-1 confined-particle transition. Recently, Skolnick ct a/. 18 confirmed that a peak 
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below the n=l electron-heavy-hole transition is phonon replica by observing the cou¬ 
pling at both the GaAs- and InAs-like bulk LO-phonon energies of InGaAs QW. Our 
measurements support these arguments. 

2. Time-resolved 


Time-resolved PL intensities of the MQW structure at 4.3K were measured over 
the spectral range from 720 lo 780nm. Lour representative streaks at given selected 
wavelengths are shown in Fig.4. The emission centered at 770nm arises from the 
recombination of band edge electron and heavy-hole. The radiation with wavelengths 
less than 770nm are from the recombination of hot carriers. The time-resolved 
luminescence intensity of the phonon replica is centered at 780nm. Each lumines¬ 
cence profile is an average from 10 individual shots using a prepulse 19 for averaging. 
The left peak of the dotted curve is the prepulse which reflects the 2ps temporal reso¬ 
lution. The right peak on the same curve is the Rayleigh scattering light from the 
sample surface which defines the "zero" time for our analysis. Several features 
appear in the data displayed in the Fig.4. (i) The rise time for all the luminescence 
profiles are instrumental; a 2ps up-limit reflects the rapid thermalization, capturing, 
and initial cooling, (ii) The shape of rise-part of luminescence profile at 780nm is 
similar to that at 770nm, but delayed by ~ jps. This implies that the emission at 
780nm does not originate from the same band as the emission at 770nm; otherwise, 
the rise-part of the luminescence at 780 should start at the same zero point as that at 
770nm. The delay of '''dps is consistent with our assignment of the emission band 
center at 780nm to the stimulated phonon-replica. This dps is the time required to 
establish a large population of NE phonons at q r . (iii) The luminescence decay time 20 
of 30ps (broken curve) for the phonon replica emission at 780nm is shorter than the 
luminescence decay time of 60ps at 770nm. This supports that the emission at 
780nm can not be due to impurity emission at high excitation. 











To study the NE phonon effect on energy relaxation process quantitatively, car¬ 
rier temperature and density as function of time should be determined simultane¬ 
ously. An expression in the time domain is introduced to fit the time-resolved PL 
data by using two adjustable parameters, namely, the carrier density 21 (n c ) and the 
electron temperature ( T c ). For the direct optical transition in MQW, the lumines¬ 
cence intensity is given by: 

I(E,,t ) = CM-e~" r ')[\M e - hh I 2 p, / , p hh f H + I 2 P e f e P,„f h Ml) 

where the p rjlh ji, -~m cjlhjh are the density of state for electron, heavy-hole, and 
light-hole; the A/ r _ /lA (//l ) is the matrix element 22 for electron to hhUh) transition; C, 
absorbed all the constant factors including the corrections for detector response and 
the transmission of each narrow band filters used; r, is the rise time of luminescence 
and is set to be ~ lps which is an up-limit of thermalizalion time of electron-hole 
system; 



is the Fermi-Dirac distributions for the electrons in conduction band (.with subscript 
e) and for the holes in valence band (with subscript h). 

An unique set parameters of T c (t ) and u, (r ) was determined by consistently 
fitting all the luminescence profiles detected at different photon energies E ,. Three 
calculated luminescence profiles corresponding the experimental data are shown in 
Fig.4 by the thick solid curves. 

The determined T c as function of time is plotted as the solid curve in Fig.5. The 
shaded area reflects the extent of the uncertainty in deducing the carrier temperature 
within the first 4ps due to our limited time resolution. The plotted data in Fig.5 is 
interpreted as follows. Carrier-carrier collisions quickly leads to a thermalized distri¬ 
bution at very high T c within lps due to a large number of hot carriers excited. 23 
Although an initial thermalization process of carriers can only be probed by fern- 
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tosecond spectroscopy. An initial cooling of the thermalixed distribution is stu¬ 
died with our present time resolution of 2ps, providing information about 
the energy relaxation and the density decay of photogenerated carriers. The 
initial cooling within first 5ps is very fast (250K/ps) reflecting both the screening 1 
and NE phonon effects are small. After the initial cooling, a large number of longitu¬ 
dinal optical phonons in a finite wavevector space accumulates due lo the finite pho¬ 
non lifetime. 24 The phonons emitted by the hot electrons are reabsorbed by the elec¬ 
trons as a reverse process of the emission giving rise to a slower cooling for the hot 
carriers. 3 It should be emphasized here that the rapid initial cooling rules out the 
importance of significant screening of electron-phonon interaction in the present 
study. If the initial screening was important, 1 the initial rapid cooling would be 
slowed and a NE phonon population would not be built up. 

Another important observation displayed in the Eig.5 is that the time constant 
for the slow decay component of carrier temperature cooling curve is ~30ps which is 
same as the decay time of the NE phonon stimulated phonon replica emission at 
780nm (see Eig.4). This suggests that ~5ps is required for the electron sys’em and 
the NE phonon system to be essentially equilibrated with each other. The coupled 
electron-phonon system decays with a common decay constant of about 30ps. It 
should be emphasized that this time constant (30ps) is the lifetime of NE phonons 
which differ from the equilibrium phonon lifetime of ~7ps. The longer lifetime for 
the NE phonons than equilibrium phonon is due to the coupling between the hot elec¬ 
trons and NE phonons. 

The determined quasi-Eermi-encrgies for electrons ar ' holes are plotted in Fig.6. 
The changes of u e and u h are very rapid with the first lOps reflecting a rapid decrease 
of the carrier density. The behavior of the degeneracy for electrons and holes with 
function of time is reversed due the difference in effective masses for the electrons 


and holes. 
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The quasi-equilibrium distribution functions for the electrons and the holes can 
be experimentally determined from T c (t ) and u r M (t ). Using these distributions the 
energy loss rates for the electrons and holes can be obtained. The energy loss rate 
(P r j, h ) is defined as follows: 


Pc . hh — 


d <E >, 


d r /„ € f {T c Me M ,e)de 1 
dt r " ■ : J 


(3) 


dt dt f, f (T c ,u, M ,e)de 

When the electron-hole system is treated as the Maxwcll-Boltzmann gas. the energy 

loss rate for the carriers (c=e,hh) is given by: 


P, = K 


. dT c it ) 


dt 


(4) 


It should be pointed that the energy loss rates for the electrons and heavy-hole 
obtained by Hq.(3) are the net loss rates. These include the energy loss due to 
electron(hole)-phonon interactions and the energy exchange between electrons and 
holes by carrier-carrier scattering. The hole-phonon scattering rate was calculated 25 
and measured 6 to be 2.5 to 3 times larger than the rate for the electrons due to the 
additional coupling through the deformation-potential. The energy transfer from 
electrons to holes should be expected in our photogencrated carrier system as long as 
the characteristic time for the electron-hole scattering is shorter than the electron- 
phonon scattering time. 


The determined energy loss rates using Eqs. 3 or 4 are shown in Fig.7. The solid 
and dotted curves were obtained from Eq.3 for the electrons and heavy-holes, respec¬ 
tively. The dot-dashed curve was obtained from Eq.(4). It should be pointed out that 
the curves obtained from either fiqsJ or 4 are the experimental data because the dis¬ 
tribution function f {T c ,u r/ ,) was experimentally determined. The broken line in 
the Fig.7 is calculated based on the simple theory 5 where the lattice is treated as a 
heat bath for the quasi-equilibrium carriers. 

The salient features of the data displayed in the Fig.7 are: 






- 15 - 


(i) The energy loss rale for the electrons is more than two orders smaller than 
predicted by the simple theory 5 (broken line). The reason for the difference arises 
from the re-heating of electrons by the NE phonons. The energy loss rate determined 
for the electrons by scattering with LO-phonons via the wave-vector dependent 
Frohlich interaction may be even smaller than what we have determined because the 
net energy loss rate determined from Eq.3 is sum of rates due to electron-phonon 
interactions and electron-hole interactions. The smaller P hh than P c is because of the 
presence of light-hole population and the different masses for the electron and hh 
which leads to different behavior of quasi-Fermi energies for the electron and hh as 
function of time. 

(ii) The net energy loss rates for the electrons and heavy-holes almost follow the 
rate determined from Eq.(4). The reasons for this are the rapid decrease of u r as 
shown in the Fig.6 as well as the slow cooling of the carrier temperature which 
result in u c / KT C <0.5, making the electron-hole system behave more like a 
Boltzmann gas. 

(iii) The energy loss rates increase as T c increases. The curves in the Fig.7 even 
bends over and approaches the maximum value of 5.3xlO -8 W for bulk GaAs. This 
"bending" behavior reflects the initial rapid cooling of hot carriers and was also 
observed by Shah et aL 6 This may be the indication for the presence of NE phonons in 
the system. 

The average phonon emission time 7,8 for hot carrier to emit a LO phonon is 
defined by: 

r„„ m ^_. (5) 

The value of r ol , ? as function of the carrier temperature is plotted in Fig.8. This is 
just another way to describe the energy relaxation process. As can be seen the r ayg is 
not a constant but carrier temperature dependent and hence time dependent. 
For T c >\200K or t < 5ps the r ayg remains at a value of about lps. At the initial 
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stage of carrier cooling only an up-limit of 0.9ps was determined from the data 
shown in the Fig.4. The theoretical predicted value should be ~-0.16/>5 7 which is 
consistent with our value within the limited time resolution. As time t > 5ps (or 
T c < 1200A: ) the T n ,. g increases quickly with decrease of the carrier temperature 
reflecting that a large NE phonon population is built up and the carrier cooling is 
suppressed. At T c = AWK , r„ l? is lOps. Ryan ct a/. 8 obtained a constant value of 
7ps for 7a\g by studies of time-resolved PL of modulation doped-MQW. The reason 
for this is simply because their time resolution of ~20ps did not allow them to mon¬ 
itor the entire carrier cooling process. 

Information about the time dependence of initial carrier population can be 
obtained using the experimentally determined distribution function f fJt . The densi¬ 
ties are obtained by the following expressions: 

n c (t ) = f it p, f c d e ; (b) 

n M U ) = f ft Phh f h d e i 17) 

(t ) = f^Plh f » de ' (8) 

where A£ is the energy separation between the hh- and lh-subbands at the zone 
center. The carrier densities for the electrons, heavy-holes, and light-holes as function 
of time are plotted in Fig.9. The quantity of In [n,. (r )] is also plotted in this figure as 
a broken curve. 

The salient feature of the data in the Fig.9 is that the carrier density decreases 
non-exponentially and varies (by a factor of 10) rapidly within the first 30ps after 
the excitation by laser pulse. One can not define a single carrier lifetime due to the 
non-exponential nature of the density curves. However, an effective carrier depletion 
time [density decreases by a factor of e~‘ from n c (t =0)] is deduced to be as short as 
lOps. The slow component of density decay curve reflects part of bimolecular recom¬ 
bination process. The fast decrease of carrier density occurring in such a short time 
period cannot be accounted for by the usual bimolecular recombination, which is a 
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much slower process, on a nanosecond time scale in a bulk GaAs and about 35Ops in 
quantum well structure with well thickness of ~5nm. 26 Several other possibilities 
should be pointed out and discussed: 

(i) Rapid carrier diffusion process may be utilized 27 to shorten the lifetime of the 
carriers within an excited spot. Assuming the velocity of photoexcited carriers mov¬ 
ing in the directions perpendicular to the z-axis is on an order of 10 7 cm Is, an 
extremely small radius of about 1 um for the excitation spot would be required to 
decrease carrier density by a factor 10 in the excited region. Since the radius of the 
excitation spot was abcut 160um in our experiment, the carrier diffusion out of the 
excited region can not result in the ultrafast depletion time of lOps. 

(ii) Auger recombination resulting from strong electron-electron Coulomb 
interaction may be very efficient and causes a substantial shortening of carrier life¬ 
time at high carrier density in a quasi-two-dimensional carrier system. But from a 
comparative study, 28 Auger recombination in semiconductor MQW structure and 
bulk were found to be close to each other. Using a value of 1.5xl0 _29 cm b r for the 
nonradiative Auger rate in bulk GaAs at room temperature. 29 Auger recombination 
lifetime at carrier density of 10 19 cm" 3 should be b70ps which is about 100 times larger 
than the carrier depletion time. However, a detailed study of Auger recombination in 
semiconductor MWQ structures under the condition of high carrier density, high car¬ 
rier temperature, low lattice temperature, and presence of a large population of Mi 
LO phonons does not exist in the literature. 

(iii) Finally, let us consider what role the NF phonons can play in reducing the 
initial carrier depletion time to lOps. As discussed in the steady-state section, there is 
a stimulated phonon emission process accompanied by an intense stimulated phonon 
replica which is mainly radiated at the sample edge and propagates parallel to the 
well plane. We found that the stimulated phonon replica detected along the y- 
direction also had a decay time constant about 30ps which is same as that detected 
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along the 7-direction. Therefore, it is most likely that the initial rapid carrier deple¬ 
tion is due to the participation of the NT LO phonons which stimulates the emission 
of phonon replica. It should be noted that this stimulated process is remarkablely 
different from the regular stimulated emission process in which phonons do not parti¬ 
cipate. As in semiconductor lasers, the stimulated emission is less effective at room 
temperature. Moreover, a longitudinal optical phonon (LO)-phonon replica plays no 
significant role in the laser operation of bulk GaAs. which is expected because of the 
weak electron-phonon coupling. Therefore, the prominence of phonon-assistant 
recombination in our quantum well structure is due to the presence of a large 
number of NH LO-phonons produced in the relaxation process of hot carriers. By 
employing the NE-phonon model we able to explain not only the slow carrier cooling 
hut also the ultrashort carrier depletion time m highly photoexcited undoped quantum 
well structures. 

V Conclusion 

With 2ps time resolution we are able to study the initial energy relaxation of 
the hot carriers and the decrease of carrier density simultaneously. The existence of a 
large population of NI: phonons in highly excited semiconductor quantum well struc¬ 
tures is experimentally verified. Its effect on the energy relaxation is to slow down 
the cooling rate after an initial rapid cooling (0-5ps). A new mechanism used to 
explain the ultrashort carrier depletion time — 10/>s deduced from the fitting of 
time-resolved PL profiles at different, emitted photon energies is proposed to be associ¬ 
ated with NE phonon stimulated phonon replica emission rather than other nonradia- 
tive processes. 
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Figure Caption 

Fig. 1 Time-integrated luminescence spectra at various lattice temperatures. The peaks 
A, U, and C are explained in the text. 

F : ig.2Triangles are the ratios of the peak intensity of B and A indicated in the Fig. 1. 
The solid curve is explained in the text. 

Fig.JTiine-integrated PI. spectra taken in the geometry as shown in the inset at 100K. 

Fig.TTime-resolved PL profiles (thin solid curves) from the symmetric MQW at 4K at 
various wavelengths. The narrow peaks within 0-1 Ops are the pre-pulses used 
for averaging. The next peak of dotted curve is the Rayleigh scattering light 
from the sample surface. The thick broken curve is generated by an expression 
of / (£ j. )=/ Jl— e ' 1T ' ]e '' Tj and the thick solid curves are generated by 
Eq.( 1). 

Fig.5Experimentally determined carrier temperature as function of time. The shaded 
area indicates the extent of uncertainty in deducing carrier temperature within 
first 4ps. 

Fig.6Experimental determined quasi-Fermi energies for electrons (up solid curve) and 
for holes (down solid curve). 

Fig.7Experimentally determined net energy loss rates for electrons (e) and heavy- 
holes (hh) obtained from two quasi-equilibrium Fermi-Dirac distribution func¬ 
tions. The dot-dashed curve is obtained from Maxwell-Boltzmann distribution 
function. The broken line is theoretical result for bulk (iaAs without the effects 
of screening and nonequilibrium phonons. The power unit is KT 1 " watts. 

Fig.BExpcrirnentally determined average phonon emission time as function of carrier 
temperature. 

Fig.9Experimentally determined carrier densities as function of time, e for electron, 
hh for heavy-hole, and lh for light-hole. The broken curve is the nature loga- 



















WAVELENGTH (nm) 









740 795 850 

WAVELENGTH (nm) 


Fig. 3 




















































Observation of Ultrafast Lateral Expansion of Energetic 


Carriers in An Asymmetric GaAs Quantum Well 

Kai Shum. H. S. Chao, M. S. Junnarkar, and R. R. Alfano 
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Departments of Electrical Engineering and Physics 
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ABSTRACT 

Lateral expansion of energetic carriers excited by an intense 0.5-ps 
laser pulse in a single asymmetric GaAs quantum well is investigated 
at 4.3K using a streak can. .ra. On a picosecond time scale, it is 
observed that the carriers expand from the excitation spot at an 
ultrafast speed of ~ \0*cm / s. A comparison, under identical condi¬ 
tions, of the asymmetric and symmetric quantum well results sheds 
light on the mechanism that leads to such an extraordinary large car¬ 
rier expansion velocity. 

PACS numbers: 66.30.Dn, 72.20. Jv., 73.40. Lq. 
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Over several years, spatial expansion of photoexcited carriers in semiconductor 
bulk and microstructures has been a subject of great interest Based on the photo- 
luminescence (PL) data, Forchel ct al. x have suggested that the fast ambipolar 
diffusion of optically excited carriers in a wide varity of semiconductors (e.g., Si, Ge, 
GaAs, CdS, and ZeSe) at a speed comparable to Fermi velocity is essential for the 
understanding the broad line shape of the luminescence spectra. Later, Steranka and 
Wolfe 2 pointed out that the fast-diffusing model 1 overestimated carrier expansion 
velocity by an order of magnitude on nanosecond time scale in indirect-gap Si due to 
neglecting the local lattice heating effect. Subsequently, Tsen 3 confirmed the results 
of Steranka and Wolfe 2 on a picosecond time scale by time-resolved Raman scattering. 
Using excitc-probe technique, Majumder cl a/. 4 found the carrier expansion velocity 
was smaller than Fermi velocity in the direct-gap CdS and CdSe. In contrast, many 
groups deduce from their experiments carrier expansion velocity v c ^ \Qpcm / s in 
GaAs 5-6 and CdSe. 7,8 This value is roughly equal to the Fermi velocity of carriers and 
exceeds group velocities of longitudinal and transverse acoustic phonons by about 2 
orders of magnitude, thus, supporting the fast-diffusing model of Forchel ct al) in 
these direct-gap polar materials. Theoretically, Wautelet and Van Vechten 9 
have inferred that hot carriers generated by nanosecond laser pulses will be 
confined near the excitation spot; while a spatial expansion of photoexcited car¬ 
riers is expected from both hydrodynamic 5 ’ 10 and thermodiffusive 11 points of 
view. Recently, Steranka and Wolfe 2 have suggested that nonequilibrium 
(NE) phonon effects 12 (phonon wind) must be taken into account in a more 
realistic description of the dynamic carrier expansion process. The NE phonons 
created in hot-carrier relaxation process provide a large driving force for car¬ 
rier expansion. 

In GaAs/AlGaAs system, the band-discontinuities at the interfaces of GaAs and 
AlGaAs leads to the formation of two-dimensional electronic energy subbands in the 
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GaAs region. Recent Raman-scattering studies by Sood ct al. 3 provided evidence for 
the presence of interface phonon modes and confined longitudinal optical (LO) and 
transverse optical (TO) phonons in GaAs well layer. One would like to kown how 
carrier confinement, modified phonon modes, and well shape influence the dynamic 
process expansion of photoexcited carriers in a quantum well structure. Most 
recently, Tsen and Morkoc 14 have shown no influence of the effects of reduced dimen¬ 
sionality for carriers and phonons on photogenerated carrier expansion in 
GaAs/AlGaAs symmetric MQW structures using picosecond time-resolved Raman 
scattering at a low carrier density of ~-10 !(j cm _i . Their measured expansion velocity 
of 8.6xl0 6 cm / s is consistent with the fast-diffusing model of Forchel ct alJ 

In this letter, we report on the first real time investigation of ultrafast hot car¬ 
rier expansion process associated with an ultrashort carrier depletion time of ~ 
lOps in an asymmetric GaAs quantum well at 4.TK using a 2-ps-resolution streak 
camera system. By studies of time-resolved PL arising from the recombination of 
energetic electron and hole detected from a lateral displaced region ~~ 3t K)um away 
from the center of the excitation spot which is ~ 160 um in diameter, it is observed 
that the hot carriers expand from the excitation spot at an ultrafast speed of ~ 
10 8 cm / s. Through a comparison of the symmetric and asymmetric quantum 
well results under identical conditions, it is concluded that the observation of such a 
fast hot carrier expansion process is made possible only in the asymmetric well 
with the higher energy barrier such as to prevent from carrier expansion in the well 
direction. Density and cooling of photogenerated carriers are independently deter¬ 
mined by analysis of time-resolved PL from the excitation spoL We have demon¬ 
strated that cooling of hot carriers is severely slowed by the NT phonon heating effect 
and the carrier density decays very rapidly within the first 30ps after the excitation 
by 0.5ps laser pulse due to the NL phonon stimulated phonon replica. 


In our experiments, an intense laser pulse of 0..“5ps at (>20nm was used to excite 
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the electron-hole pair density of ~ lO'^cm -3 . The excitation source is from a 
colliding-pulse passive-mode-locked dye laser amplified by a four-stage dye amplifier 
pumped with a frequency doubled Nd: YAG laser operated at 20Hz. The undoped 
GaAs/Al^Gaj x As asymmetric and symmetric quantum wells investigated in this 
work were grown by molecular-beam epitaxy on a (001 Loriented undoped GaAs sub¬ 
strate. The asymmetric quantum well consists of layers from the front surface to 
substrate 150A Al 03 Ga 07 As, 50A GaAs well, 200A AlAs (asymmetric barrier), and 
followed by 1 um Al 03 Ga 07 As, and 0.2 um GaAs buffer. The symmetric well struc¬ 
ture used for the comparative studies consists of 50 periods of 55A GaAs and 100A 
Al 03 Ga Q7 As layers followed by a 1.2 um GaAs buffer layer. The excellent quality of 
the samples are confirmed by well resolved heavy-hole and light-hole excitonic struc¬ 
tures in room temperature PI spectra. 

First, we discuss the measurements of the differences in the spatial widths of 
time-integrated laser and PL. intensity profiles on the sample surface. 1 he hoi PI. 
profiles (solid curves) detected within a spectral window of 35nm ( ( K)meV ) centered 
at 700nm (1.771eV) for the asymmetric and symmetric wells along with the profiles 
(broken) of corresponding exciting laser at a fixed intensity are plotted in Fig.la and 
Fig.lb, respectively. 1 he gain of detecting system was adjusted such that the max¬ 
imum intensities recorded for the PL and laser w r ere roughly same. The insets show 
the potential profiles of the asymmetric and symmetric wells. It is clear that the 
half-width of PL profile for the asymmetric well is much broader than for the sym¬ 
metric wells. This implies that the energetic carriers (about 3 LO-phonon energies 
above the subband edge) in the asymmetric well is much more mobile than in the 
symmetric well. 

By making use of the two electronic windows of the streak camera, time- 
resolved PL intensities of the asymmetric well, I w] (t) and I w2 (t), from the excited 
and unexcited regions, respectively, were measured. The subscripts wl and w2 






(r i“*r i) stand for the two windows as being positioned in the Pig.l. The l u-1 (t) and 
I w2 (t) (thin solid curves) measured at 4.TK with ~2ps time resolution are displayed 
in F ; ig.2. The peak of the broken curve is the Rayleigh scattering light (620nm) from 
the sample surface which defines the "zero" time for our analysis. The thick solid 
curves are the theoretical fits which will be discussed below. As can t>e seen the rise 
of the l u .j is time-: ^solution limited, while the rise of l w1 is much slower than that 
of the I^.j. The slow rise time 15 of I w1 (r,“ 2 ^ lOOps) is essential to suppoort that 
the luminescence detected from the laterally displaced w2 indeed arises from the 
recombination of electron and hole which come from the wl. It eliminates one possi¬ 
ble artifact that may lead to anomalously strong luminescence detected from the w2. 
Such a possibility arises from the 1 uni layer of AI (n Ga ()7 As which is next to the 
AlAs barrier. This layer is also highly excited by the (>2()nm laser pulse. Possible 
recombination radiation (although we did not detect for the asymmetric well, neither 
for the symmetric well) from Al ()3 Ga 07 As layer will excite carriers within the well. 
This secondary excitation would not be limited to the originally excited region giving 
rise secondary luminescence. Since the speed of light in GaAs is about l(X)um /ps, the 
luminescence intensity I w2 would rise quite similarly with the I wj . This is contrary 
to our experimental observation. 

Having known the expansion distance (294 j/ m ) from the center of the w 1 to the 
inner edge of the w2 and assuming the time it takes for carriers to expand from the 
wl to w2 where they recombine is equal to the rise time of l w1 , an up-limit of the 
expansion velocity t>, is estimated to be about 2.94xlO s cm / * {29Aum / lOOps). 

The decay of the l wJ is dependent upon a transient carrier distribution in the wl 
which determines the number of carriers in the given energy space as function of 
time. In order to obtain the carrier distribution (assumed to be hermi-Dirac distribu¬ 
tion), the time-resolved PI. profiles detected in narrow (20meV) spectral regions 
which cover the whole emission band were measured from the region wl at 4.3K. 
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We have introduced an expression 16 using the carrier temperature [7\ (t )] and the 
quasi-Fermi-energy for electrons (; )] as a set of parameters to consistently fit these 
measured profiles, three typical measured profiles as well as theoretically fitted 
curves (thick solid) are plotted in T'ig.Ta. The subband edge luminescence is centered 
~760nm. The cooling of hot carriers is substantially slowed down after ~5ps due 
to reabsorption of NT phonons by the carriers. 16 ,8 The electron density n c it ) and 
In [n,. (r )] are plotted in Fig.ib. The salient feature of the plotted data is that the car¬ 
rier density decreases nonexponentially and very rapidly form 0 to 30ps. One can not 
define a single lifetime to describe the carrier density decay. However, an effective 
carrier depletion time [the time for n,. to decrease by a factor of e -1 from n,. (t =(>)] 
can lx* deduced to be as short as lOps. Such an extremely rapid carrier density 
decrease can not lx- accounted for by an usual bimolecular recombination which is in 
nanosecond lime scale for a bulk Ha As and about 3.M)ps lor .1 quantum well structure 
with well thickness of ~5()A. 1<y But it can be explained consistently b\ the participa¬ 
tion of NT phonons in the highly excited system. Since NT phonons aie locali/ed both 
in a wave-vector space and in a real physical space (in the well plane) they will 
behave like coherent Bosons. The existence of the coherent Bosons will tend to 
increase the number of the Bosons in the system with a rate proportional to the 
present number of Bosons. This is a stimulated phonon emission process and can be 
accomplished by a phonon replica, when an electron recombines with a hole, not only 
a photon but also a phonon or more phonons will be emitted to join Ifoson system. 
Such a NT phonon stimulated phonon replica has been reported previously by 
Holoayak ct aL~° We have recently confirmed 16 this novel process in the symmetric 
wells by monitoring the phonon replica which is mainly radiated in the lateral x-y 
plane in time domain as well as detailed polarization studies of integrated lumines¬ 
cence spectra. 

Using the determined carrier distribution, the solid curve as shown in Fig.2a is 
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Using the above expression to fit the measured I .two parameters of I)' and r/ 2 
have been precisely determined to be 0.7x1 (/'cm ■’/ .» and 13ps, respectively. 

It should be investigated whether the determined values for D’ and ry 2 can be 
interpreted from the t herrmxlifTusive model. 1,11 Several aspects are discussed below: 

(1) The value of O' is four orders magnitude greater than the regular impurity 
scattering limited diffusivity. The enhancement of /)' may arise from high carrier 
temperature and the screening effects 8 at high carrier density. However, the screen¬ 
ing effects do not seem to be effective since the initial carrier cooling is very rapid. 16 















(2) The decay rates of carriers which are located in the same energy space but in 
the different spatial regions of the w 1 and w2 are comparable reflecting the dominant 
process for carrier density decay is not due to the carrier expansion out from the wl 
but due to the NE phonon stimulated phonon replica. It is further inferred that a 
large number of the NT phonons emitted by hot electrons may exist in an area of at 
least ~ 300 um in diameter. These phonons which have momenta parallel to the x-y 
plane relax to acoustic phonons of large momentum and then these acoustic phonons 
may drive 2 the carriers outside the excited region. A rough estimation seems support 
this argument. Let us take a value of 10 -8 dynes per pair 2 for acoustic phonon-wind 
force, an excess energy of 140meV per electron-hole pair, and an acceleration time of 
13ps. 1'hen, a value of 2.8x1 O s cm / s for carrier expansion velocity is predicted by 
phonon-wind model. 2,12 

(3) From the determined D and t} 2 , the expansion velocity is calculated by 

v,, = -v / = 2.3x10’em / s which is consistent to the thermodiffusiw model 1,11 

V Tj - 

and the estimated up-limit. In a two-dimensional electron system quasi Termi-encrgy 
is proportional to carrier density. Therefore, the Fermi-velocity v F and hence \y is 
proportional to the square root of the carrier density. Then, our observed 
v, = 2~-3xlO s cm / s at the high carrier density of ~ 10 iq cm -3 is consistent with the 
measured v p = 8.6x10 l, cm Is by Tsen and Morkoc 14 at the low carrier density of ~ 
10 lh cm“ 3 . 

(4) Using the determined D' and rf 2 , an effective "diffusion length" is calcu¬ 
lated to be 30 um . According to this value, there would be no significant carriers in 
the region of w2 which is ~3OOum away from the excitation center. But the total 
detected photons from the wl and w2 within the first 20ps after the excitation by 
the laser pulse is comparable. This might be again interpreted by the phonon-wind 
model 2 in which the concept of conventional diffusion length may not apply since 
the collisions between carriers and phonons no longer limit how far carriers can be 
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drifted. 

(5) Finally, it must be noted that the anomalously fast expansion is only 
observed in the asymmetric well. This is very interesting because such a simple 
modification of quantum well shape in comparison with the symmetric well makes a 
big difference on the dynamic carrier expansion process. Due to the asymmetric well 
shape there are no transmission states (virtual states or unconfined states) above the 
lower energy barrier unlike the symmetric well case. This is an intrinsic property of 
the asymmetric quantum well. The photogenerated carriers may be immediately 
scattered by the higher energy barrier and obtain momenta from the interface which 
assists carriers to expand rapidly in x-y plane. 

In summary, we have observed for the first time the ultrafast expansion of the 
energetic carriers in the asymmetric quantum well. It is found that the unique well 
shape of the asymmetric quantum well is responsible lor the observation of the 
extraordinary large expansion velocity of the energetic carriers. However, no distinc¬ 
tion has been made among various models to explain the observed ultrafast carrier 
expansion process. Both the rapid carrier density decay and the slow cooling of hot 
carriers in the asymmetric and symmetric wells are attributed to the presence of the 
nonequilibrium phonons generated by the energy relaxation of photoexcited carriers. 

The research was supported by the U.S. Air Force Office of Scientific Research 
under Grant No. AFOSR-86-0031. We thank K. Bajaj for helpful discussions. 
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Figur e Ca ption 

Fig. 1 Time-integrated spatial profiles of laser and the luminescence in a spectral win¬ 
dow of 35nm centered at 700nm at 4.3K from the asymmetric well (a) and 
symmetric wells (b). The wl ad w2 are regions where the time-resolved 
luminescence profiles are obtained. The insets show the potential profiles of the 
asymmetric and symmetric wells. 

Fig.2Time-resolved photoluminescence of the asymmetric well at 4.3K from the spa¬ 
tial region of the wl (a) and w2 (b). The solid curves are the theoretical fits as 
discussed in the text. 

Fig.3(a) Time-resolved photoluminescence from the asymmetric quantum well at 
4.3K at various wavelengths, (b) The electron density decay curves. 
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Electron-hole recombination lifetimes in a quasi-zero-dimensional electron 
system in CdS^Se^, 
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The recombination lifetimes for the radial and angular quantum number conserved 1S-1S and 
1 P~\P transitions from three-dimensionally confined electrons in CdS, Se, _ , were measured 
by time-resolved photoluminescence (PL). The assignment of the observed transitions was 
supported by calculations of eigen energy levels and squared matrix element ratio for these 
transitions as well as well-resolved PL peaks arising from IS-15 and \P-\P transitions. 


Recently, carriers localized in semiconductor micro- 
structures have attracted much attention because of their 
novel optical properties for potential device applications. 
Electrons localized in semiconductor crystallites of diame¬ 
ters ranging from 30 to 800 A embedded in a transparent 
insulating matrix are confined in three dimensions. 1 ' 3 Such 
quasi-zero-dimensional (quasi-OD) electron systems also 
exist in bulk alloy semiconductors due to compositional fluc¬ 
tuations. 4 Electronic motion in these systems no longer fol¬ 
lows a well-defined energy momentum relation because the 
Hamiltonian of the system does not commute with the mo¬ 
mentum operator due to potential discontinuity at the crys¬ 
tallite surface. The envelope wave functions in spherical co¬ 
ordinates [ r,6,d 1 and eigen energies of conduction electrons 
and valence holes localized in an infinite spherical well with¬ 
in the effective mass approximation are given by 

\nLm) - C„J L (x*L r ) Lm ) (1) 

and 

E mt =frx z nL nm th a\ (2) 

respectively, where the subscripts n, L, and m are effective 
radial, angular, and magnetic quantum numbers, respective¬ 
ly; the j and Y are the spherical Bessel and spherical Har¬ 
monic functions, respectively; a is the radius of crystallite; 
and m, k is the effective mass of electron or the isotropic hole 
mass. 5 The values of x*l for the lowest two states of either 
conduction electron or valence hole IS and IP are n and 
4.49, respectively. S stands for L = 0 and P for L — 1. The 
allowed transitions which conserve angular and radial quan¬ 
tum numbers are 1S-1S, 1P-1P, and higher transitions. The 
physical picture of quasi-OD electron system described 
above has been experimentally verified by several groups. 1 ' 3 
Large and fast optical nonlinearities arising from photogen¬ 
erated electron hole pairs have been reported in quasi-OD 
electron systems. 6 7 

In this letter, we report on measurements of ultrashort 
recombination lifetimes of the radial and angular quantum 
number conserved lS-lSand \P-\P transitions in quasi-OD 
electron systems in CdS, Se, _ „ at 4,3 K using a streak cam¬ 
era detection system. The time-resolved photoluminescence 
(PL) detected at various emission energies allows us to un¬ 
ambiguously identify the IS-IS and 1P-1P transitions. The 
recombination lifetime of IP-IP transition was measured to 
be 3.5 times shorter than the IS- IS transition. The ultrafast 
decay of the IPexcitation may have practical importance for 
the construction of ultrafast reversible optical switches. The 


assignment of the observed S and P transitions was support¬ 
ed by calculations of eigen energy levels and squared matrix 
element ratio for these transitions, as well as the observation 
of optical transitions in steady-state PL. 

The samples investigated were four optical glass filters 
Coming 2-61, 2-59, 2-58, and 2-64. The value of x for each 
sample was accurately obtained from chemical analysis. 8 A 
second harmonic (530 nm) of a Nd-glass laser pulse of 8 ps 
duration was used to excite the samples on the front surface. 
The maximum optical energy incident onto the front surface 
was about 40/zJ. The spot size was about 8x 10' 3 cnr. For 
steady-state PL experiments, an argon-ion laser, a Spex dou¬ 
ble grating spectrometer, an S-20 photomultiplier, and a 
lock-in amplifier were used. 

The time-resolved PL profiles obtained at various emis¬ 
sion energies at 4.2 K for sample 2-58 are shown in Figs. 
1(a)—1 (e). The left curve in Fig. 1(a) shows the temporal 
profile of the exciting laser pulse which reflects the time reso¬ 
lution of the detection system (10 ps). The dotted curves in 
Fig. 1 are the double-exponential fits to the data with a value 
of 12 ps for the rise time r,. The decay times of these time- 
resolved luminescence profiles show only two distinct val¬ 
ues. For emission energies ranging from 2.213 to 2.175 eV 
and from 2.1 to 1.967 eV the decay times are 29 and 100 ps, 
respectively. 

The decay time as function of emitted photon energy E 
for four quasi-OD electron systems at 4.2 K is summarized in 
Fig. 2. For a comparison, the exciton lifetime versus exciton 
energy in a quasi-OD system in a bulk CdSosjSe^ alloy 
compound studied by Kash el al. 9 is also included in this 
figure. The most remarkable feature of interest in the data is 
the appearance of two distinct energy regions in which the 
decay time is different by a factor of about 3.5 for the samples 
2-64, 2-58 and the bulk CdSosjSeo,, while the other two 
samples do not show this feature because either IP electron 
and hole were not excited or the luminescence from the 1P- 
1P electron-hole recombination was not detected. It is ex¬ 
pected that the decay time should be nearly invariable over 
the energy range of emission if only the lowest confined state 
IS for electrons and holes is occupied. The wide spectral 
range most likely reflects the crystallite size distribution and 
the fluctuations in the value of x from crystallite to crystal¬ 
lite. When the two lowest states IS and IP are substantially 
occupied and the recombination lifetimes associated with 
these two states are considerably different, a steplike change 
of decay time will be observed in two distinct energy regions. 
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FIG l. Time-resolved PL profiles de¬ 
tected at different energies (a)-(e) for 
the sample 2-58 at 4 2 K. The dotted 
curves are the theoretical fit with 12 ps 
for rise time and decay times are indi¬ 
cated on the corresponding profiles. 
The left curve of (a) is the temporal 
response profile to the laser pulse used 
to excite the samples. 


The exact lifetime ratio in the two energy regions depends on 
the transition matrix elements. The above argument should 
in principle explain what we observed. However, the energy 
dependence of localized exciton lifetime in the bulk com¬ 
pound CdSo „ Se<j 47 as displayed in Fig. 2(e) also exhibits a 
steplike feature. The explanation given by Kash et al? was 
based on the model suggested by Cohen and Sturge 4 where 
the exciton migration goes from a site with higher energy 
into another site with lower energy. This raises a question 
whether the exciton migration mechanism can apply to the 
present case. Since exciton migration must involve a trans¬ 
port over a relatively large distance due to the nature of the 
exciton-phonon interaction, 4 such migration may occur in 
relatively large crystallites. This apparently may explain the 
results observed in samples 2-58 and 2-64. In smaller crystal¬ 
lites the migration does not occur in correspondence with 
the samples 2-61. However, the exciton migration picture 
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FIG. 2. Lifetimes »re plotted as function of emitted photon energy. The 
sample temperature for (a), (b), (c), (d) is 4.2 K and for (e) is 2 K. The 
arrows indicate the energy positions for the 1S-I5 transition (S’) and the 
\P-\P transition {P). 


cannot explain the steplike energy dependence oflifetime as 
well as the ratio with the value of about 3.5 because the depth 
of a potential well caused by compositional fluctuations is 
entirely random in nature. Furthermore, the exciton lifetime 
would increase monotonically with a decrease of exciton en¬ 
ergy. 

A two-state (15, IF) model is introduced here which 
consistently interprets all the results shown in Fig. 2 as well 
as other results observed by Cohen and Sturge 4 and reported 
by Kash et al? It is possible that the excitons are confined by 
the potential wells in the bulk alloy compound which have 
two lowest states designated by 15 and IF. The exciton in 
higher lying state IF can be scattered to the lower lying state 
15 by emission of phonons instead of migrating from site to 
site. 

In order to verify the two-state model, we have calculat¬ 
ed the eigen energy levels of the quasi-OD electron systems 
and the recombination lifetime ratio for 15-15 and 1F-1F 
transitions. The confinement energy for the ground states 
(£ is ) is equal to E + E * s for the quasi-OD electron sys¬ 
tems, where E' s and are the lowest confinement ener¬ 
gies for electrons and holes, respectively. These energies 
were calculated from the bulk band gaps 10 at 300 K and the 
measured peak energies of the first derivative of room-tem¬ 
perature reflectance (dR/dA) by subtraction. Using the 
electron effective mass, the isotropic hole effective mass giv¬ 
en in Ref. 4, and the measured the effective diameter of 
crystallite for each sample and eigen energies were calculat¬ 
ed based on Eq. (2) and listed in Table I. The emitted photon 
energies at 4.2 K for 15-15and 1F-1F transitions in Table I 
were obtained by adding corresponding confinement ener¬ 
gies to the values of bulk gap 10 at 4.2 K. These calculated 
energy positions are located in Fig. 2 as arrows labeled by 5 
and F for clarity and are in reasonable good agreement with 
the two-state model. 

The radiative lifetime ratio for 15-15 and IF-IF can be 
predicted following the approach given by Casey and Pan- 
ish." The matrix element for a transition from localized 
state in conduction band to the localized state in valence 
band is M = M b M„ y , where M b is the average matrix ele¬ 
ment for the Bloch states for bands in absence of new eigen 
states due to the confinement and M my is the envelope part 
of the matrix element which is M my = (1Z. |\L ). Since the 
IF state is a threefold degeneracy state with magnetic quan¬ 
tum number of — 1,0, and 1, the envelope matrix element 
ratio of |M tnv |^_ 1P to |A/ CT> |} s _ 1S can be readily calculat¬ 
ed to be 5. When the 4 K recombination lifetimes are purely 
radiative, then the ratio for r, s /r, P should be equal to 5. The 
discrepancy between the predicted value of 5 and the experi¬ 
mentally measured value of 3.5 for r 15 /r tP arises from the 
existence of nonradiative transitions. A simple rate equation 
analysis helps to substantiate this statement. Suppose 1/ 
r„ = 1 /t„ + 1 /r„„ where i = 5,F; r n(p) and r vlp) are the 
radiative and nonradiative recombination lifetimes for the 
15(F)-15(F) transition, respectively. The ratio is 
r is/Tif = (r n /T rp )[(T rp +T„ p )/(T n +r„)]. When 
~r np and •< r„ .r,,,, then r ls /r IP -1; hence, there should be 
only one lifetime. This is opposite to the experimental re¬ 
sults. On the other hand, the ratio T xs /r ir is about 3.5 when 
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TABLE I Measured energy peak of first derivative of reflection at 300 K. bulk energy gaps at 4 2 and 300 K, the calculated confinement energies, emission 
energies for 1S-1S and \P-\P transitions and effective diameters (d). the measured lifetimes for I S-15and \P-\P transitions as well as their ratio are 
displayed 
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‘ From Ref 10 

"Measured energy separation (meV) between 5 and P peaks at 300 K from Fig. 3(a). 
‘ Peaks are broad and unresolved. 

d The anticipated values assuming a same ratio t u /t ip for the samples 2-58 and 2-64 


nonradiative recombination lifetimes (r„, = r ) are about 
5 times larger than the radiative lifetime r rp of the excited 
state IF for the assumption of r [P ~r and r is ~r„. 

Steady-state PL studies at 4 and 300 K were performed 
to substantiate our findings from the picosecond studies. In 
Fig. 3(a), 300 K PL spectra of sample 2-61 at various excita¬ 
tion levels (see figure caption) are plotted. The broken verti¬ 
cal line indicates the peak position of dR /dA. As expected 
the salient feature in the luminescence spectra is the appear¬ 
ance of a two-peak structure separated by 220 meV on the 
high-energy side of the broken vertical line. Based on the 
two-state model we attribute these two peaks to the IS-15 
and 1 P- 1 P optical transitions’ as indicated by S and P in the 
figure. 

The 4 K PL spectrum of 2-61 excited by the 457.9-nm 
line of the argon-ion laser is depicted in Fig. 3(b). As can be 
seen from the spectrum, the position of the 15-15 peak shifts 
to higher energy and its FWHM (58 meV) decreases by a 


51; 

(0)2-61 

300K 


i 
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WAVELENGTH(nm) 

FIG. 3. (a) Photoluminescence spectra of sample 2-61 at various excitation 
levels taken at room temperature using the 488 nm line of an argon-ion 
laser. The excitation levels (intensity scale) were 2.5 (2) and 10(10) times 
higher for the broken curve and the dot-dashed curve than the solid curve, 
respectively. The broken vertical line indicates the peak position of the first 
derivative of reflectance at room temperature 5 and P stand for the 15-15 
and IP-IP transitions, respectively (b) Photoluminescence spectra of2-6l 
at 4.3 K. The inset shows that the 15 peak broadening is dominated by the 
diameter fluctuation of microcrystallites. The solid dots are calculated by 
assuming a Gaussian shape of the probability density function for the ran¬ 
dom variable d with a variance of 25 A : . 
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factor of 2 in comparison with room-temperature data. The 
latter indicates that the broadening of the 15-15 transition at 
room temperature is not entirely dominated by the size fluc¬ 
tuation. The inset of Fig. 3(b) shows a spectrum fit by as¬ 
suming a Gaussian shape of the probability density function 
for the random variable d. The fit gives a mean value of 74 A 
for the effective diameter and a full width at half-maximum 
of 15 A which characterizes the size fluctuation of.micro¬ 
crystallites. The structure next to the main 15 peak is from 
the IP-IP transition since its energy separation with the IF 
peak is about 210 meV. 

We have also measured 4 and 300 K luminescence spec¬ 
tra of 2-59, 2-58, and 2-64 samples. The two-peak structure 
was observed for 2-59 at both 4 and 300 K. The peaks were 
broader for 2-58 and 2-64 than for the other two samples and 
the IF structure was not clearly identified. 

In summary, we have reported on the measurements of 
recombination lifetimes of 15-15 and 1F-1F transitions in 
quasi-OD electron systems in CdS.Se, _ *. The assignment 
to these observed transitions is supported by calculations of 
eigen energy levels and the corresponding squared matrix 
element ratio as well as the direct observation of optical tran¬ 
sitions between quantized energy levels 15 and IF in the 
conduction band and the valence band. 

This work was supported by Air Force Office of Scien¬ 
tific Research AFOSR-86-0031. 
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Asymmetric GaAs Quantum Well3 
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Departments of Physics and Electrical Engineering 
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Subpicosecond time-resolved pnotoluminescence studies of a symmetric multiple quantum well structure and 
an asymmetric quantum well with a similar well thickness have revealed several important experimental 
results on the behavior of photoexcited carriers in these microstructures: (1) The energy relaxation process 
is substantially suppressed cue to the existence of a large population of ncnequilibrium phonons after an 
initial rapid cooling. (2) The two different well structures do not show any difference in the energy 
relaxation process under the same experimental conditions. (3) The photoexoitec carrier density deduced 
from exact and consistent fittings of the time-resolved pnotoluminescence profiles at various emitted pnoton 
energies decreases nonexpcr.ent; al iy and very rapidly within the first 30 ps after the end of a 0.5 ps pulse. 
An effective carrier depletion time is determined to be as short as 10 ps in either well. A mechanism which 
leads to such short carrier depletion time is found to be associated with r he nonequilibrium phonon enhanced 
p no non replica emission. Unlike conventional stimulated emission s<>r.av;cr, this mechanism is even more 
effective in exnaustir.g carriers from the system at higher temperatures. (a) The difference in potential 
well profile corresponds to a very different behavior cf carrier diffusion process in the lateral well plane. 
The carrier diffusion ;s enhanced m the well plane of the asymmetric well by restricting carrier diffusion 
in the growth direction. Diffusivity 0 of the pnotoexcit°J carriers :r. the asymmetric well has been directly 
determined to be 12‘cs*/s at u.JK, which :s about four cr :»rs of magnitude l arger than the value in bulk 
GaAs. 


Over the last few years, there r.as been a growing interest to understand the fundamental interactions of 
confined and free electrons and holes with ‘he phonons in bulk semiconductor and its multiple quantum well 


structures. This understanding is crucial for the operation of uitr.anigh spee; 
principles of ultrafast carrier dynamics. The slowing of carrier ceding in i 
attributed to the screening 1 of the electron-phonon interaction or the re-heatu 
phonons.' 1 ' 1 The first com pan ten* of tr.e hot-"lent nor. coding rates in un.dsped mul 1 
structures and in the bulk GaAs was achieved by time-resolv«d measurements of opti 
The cooling rates were approximately the same for bulk and MQW at a carrier densif 
is expected from a simple theory* where the lattice is ‘re ate a as a heat bath f 
carriers. Quasi-steady state experiments have teen carried out by two group 


1 devices based on the 
bulk semiconductors was 
r.g of carriers by het- 
tiple quantum well (MQW) 
.cal absorption and gain, 
y of 2.5xl0‘ 7 cm -1 . This 
'or the quasi-equilibrium 
is 6 ' 7 to investigate the 


interaction of electrons and holes with phonons in modulation-doped MQW structures. These experiments have 
generated conflicting results regarding the presence of nonequi 1 ibrium phonons. The first measurements 
of time-resolved photoluminescence from modulation dop'd semiconductor MQW with greater than 20 ps time 
resolution were reported by Ryan rt al.“ who found that the cooling of hot carriers was anomalously slow 
after R0 ps. By reviewing these and other previous works we noticed three important points: First, by 

studying the carrier dynamics in modulation-doped MQW structures which are designed for high mobility 
devices, one can obtain the energy loss rates of the photo-injected majority carriers by extracting the 
carrier temperature as a function of time from the high energy tails of time-resolved pnotoluminescence 
spectra. Due to a large extrinsic carrier density, at least comparable to the priotoexcited carrier density, 
one cannot obtain information on photo-injected carrier lifetimes. In these experiments the carrier lifetime 
is essentially close to the value in bulk which is about 1 ns.“ It 13 necessary to study photo-injected 
car-ier dynamics in undoped MQW structures in order to obtain information on both the energy loss rates and 

carrier lifetimes. The latter is cruel.1 for applications ir: high speed optoelectronic devices. Several 

attempts have been made to reduce the carrier lifetime through nonradiativ" recombination processes.’■ 10 
Second, one cannot extract information on the initial carrier relaxation process from both quasi-steady state 
experiments and time-resolved measurements with time resolution greater than ID ps. Study of the initial 
carrier relaxation process is crucial to substantiate the existence of nonequilibrium phonons and Its 

importance to the hot-carripr relaxation since the phonon lifetime is about 5-7 ps. Third, the rapid 

expansion of photoexcited electron-hole plasma at high density in quantum well structures has not been 
experimentally investigated on a picosecond time scale. Only an indirect measurement of the expansion of 
electron-hole plasma in G 3 A 3 MQW structures at low carrier density has been reported where the drift 
velocity of the electrons was comparable to the Fermi-ve 1 oci t y, ind therefore, the transport of the 
photoexcited carriers In GaAs quantum wells was thermodiffuoive. 11 

In this paper, several experimental results are reported from the measurements of time-resolved 
photoluminescence with 2 ps^time resolution as well as t:studies of lattice temperature dependence, 
excitation dependence, and polarization of integrated luminescence spectra from an undoped GaAs symmetric 
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MQW structure and an asymmetric Quantum well. We have found: (1) The energy relaxation process is 

substantially suppressed due to the existence of a large population of nonequilibrium phonons after an 
initial rapid cooling. (2) The two different well structures do not show any difference in the energy 

relaxation process under the same experimental conditions. (3) The pnotoexcited carrier density 'tracted 
from fittings of time-resolved phctoluminescence profiles at various emitted pnonon energies decreases r.on- 
exponentially and very rapidly within the first 30 ps after the <-nd 0 f a 0.5 ps laser pulse. An effective 

carrier depletion time is determined to be as she, t as !0 ps in either type of wells. A mechanism wnich 

leads to such a short carrier depletion time is associated with nonequilibrium phonon ennanced phonon- 
replica-emission. Unlike the conventional stimulated emission process this mechanism is even more effective 
in exhausting carriers from the system at higher temperatures. The bandgap renormalization in the 
photoexcited region is not observed from integrated luminescence spectra as a consequence of ultrashort 
carrier depletion time and hot carrier temperature. (A) There is no or little carrier diffusion in the 
lateral well plane: of the symmetric wells except for the carriers at the subband edge, while in the 
asymmetric well, the diffusion process is enhanced in the lateral well plane by restricting carrier diffusion 
in the growth direction. The diffusivity, D, of photoexcited carriers in the asymmetric well has been 
determined to be 10 e cm 2 /s at A.3K, which is about i'cur crcers of magnitude larger than corresponding value in 
bulk GaAs. 


Sam pies 


The undoped GaAs/Al x Gai- x As symmetric MQW structure and the asymmetric quantum well investigated in this 
work were grown by molecular-beam epitaxy on a (001 )-oriented undoped GaAs substrate. The symmetric MQW 
consists of 50 periods of 55A thick GaAs and lQoX thick Al 0 _,Ca 0 . ,As layers followed by a 1.2 urn buffer GaAs 
layer. The asymmetric quantum well consists of layers frem the front surface to substrate 150S. Al„ ,Ga 0 .-,As, 
50A GaAs well, 200A AlAs (asymmetric barrier), followed by 1 urn Al<,_ ,Ga 0 ,As, and 0.2 urn GaAs buffer" layer on 
Si-doped GaAs substrate. The lateral sizes are 2x4mm, 2x6mm for the symmetric MQW structure and tne 
asymmetric well, respectively. The samples were mounted on a cold finger in a helium cryostat. 



S' 

wS 


Experimental 


An ultrashort light pulse of 500 fs duration at 620 r.m was used to excite the electron-hole pairs with a 
carrier density of 10 15 cm -3 . This light source was generated from a collidir.g-pulse passive 1 ' mode-locked dye 
laser and amplified by a four-stage dye amplifier pumped with a frequency doubled Nd: TAG laser at 20 Hz. 
The photoluminescence was spectrally filtered using various narrow bana filters, ana temporally dispersed by 
a 2 ps time resolution Hamamatsu streak camera system. The luminescence intensities were corrected for 
nonlinearity of the streak rate, the spectral response, and the transmission of each narrow band filter. The 
time resolution for the present work is about 2 ps. The experimental apparatus used in this research has 
been described in detail by Ho et al. 12 


I 


Schematic energy band diagrams are shown in Fig. 1 for a symmetric and an asymmetric quantum well. Using 
a finite-depth square well model and material parameters given in Ref. 13 the subband energy positions in the 
wells wrre calculated. There are two subbands for electrons (?) and heavy-holes (hh) and only one subband 
for light-holes (lh). The good quality of the samples is confirmed by spectrally well resolved heavy-hole 
and light-hole excitonic structures at room temperature. The FWHMs for the heavy-hole excitons are about 15 
meV and 12.7 meV for the symmetric and asymmetric quantum wells, respectively. 


The electron-hole pairs were generated by photo-excitation. These photogenerated carriers therm 3 lize 
within themselves via carrier-carrier interaction to roach a quasi-equilibrium state in less than 1 ps. 
Simultaneously, three other processes take place. The first one is the loss of energy of the electron-hole 
system through the electron-phonon interaction. The second process is the carrier diffusion process due to 
the gradient of carrier density between photoexcited and non-excited areas. The third process is 
recombination of electron-hole pair through either nonradiative or radiative mechanisms. The work reported 
here on the photoluminescence experiments uses the last process to probe the first two processes. 





a) Steady State Luminescence 


Energy Relaxation 


V. 



1 


Time-integrated luminescence spectra of the symmetric MQW structure excited by the subpicosecond laser 
pulse at various lattice temperatures (T^) are shown in Fig. 2. These spectra were taken in the z(y,y)z 
configuratlon. The conventional Raman notation (x,y,z are the crystal axes, and z is the growth direction) is 

used. In order to eliminate radiation from the edges of the sample the luminescence spot wa 3 imaged on an 

aperture with the aid of a 3treak camera and then focused onto a vertical slit of a grating spectrometer. 
Several features are displayed in the data. 

(1) The emission peak (A) on the high energy side of the spectra is from the recrmbinatlon of 

photogenerated n-1 electrons and n-1 heavy-holes. This peak shirts towards the low energy side as T L 

Increases. A high energy tail on the main peak (A) develops with the increase of T L . For T L >150K, a 
shoulder on the high energy side of the peak A arising from the recombination of n-1 electrons and n-1 light- 
holes can be identified. 
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(2) A broad emission band (C) at the low energy side of the spectra arises from the 1.2 um GaAs buffer 
layer. The total emission intensity decreases as increases. 

(3) The most interesting feature of the data is the appearance of an emission band (3) below the n-1 
electron to n-1 heavy-hole transition which is attributed to the emission from the nonequilibrium phonon 
enhanced phonon replica. 

The following reasons support this assignment: 

(i) The B emission band does not appear in the steady-state photoluminescence spectra taken at low 
power excitation of about 1W/cm 2 at ii. 3K using a 488 nm line of an argon-ion laser. However, the 
photoluminescence spectra using a weak train of laser pulses (- 200 fs) directly from the colliding-pulse 
passively mode-locked (CPM) dye laser oscillator as the excitation source, with an excitation power density 
in the range of 10" 5 -10" J W/cm 2 , show a weak electron-acceptor emission band separated by 17 meV from the n-1 
electron-heavy-hole transition accompanied by its phonon-replica 1 ' at low temperatures. This extrinsic 
emission band disappeared completely when the sample temperature was raised to - 80K, whereas the 3 band in 
Fig. 2 exists up to room temperature. Moreover, the relative time-integrated intensity ratio between B and A 
bands in Fig. 2 detected under amplified 0.5 ps light pulse excitation decreases as the excitation power 
density decreases and the B band is not visible when the excitation power density is lower than 10"’P m , where 
P m - 10 12 W/cm 2 is the maximum value of the excitation power density. This contradicts what is generally 
expected for an impurity emission. Since the concentration of acceptors is low in our sample as confirmed by 
the luminescence studies using the weak train of 200 fs light pulses, the impurity emission should be readily 
apparent and more pronounced at lower excitation. Therefore, the B emission band shown in Fig. 2 cannot be 
attributed to electron to acceptor luminescence. 

(ii) The ratio of emission intensities of B and A emission bands increases as T^ increases. The values 
of the ratio as a function of Tl are shown in Fig. 3 as triangles. This behavior is remarkably different 
from the conventional stimulated emission process in which phonons do not participate. As in semiconductor 
lasers, the stimulated emission is less effective when the sample is at room temperature. To convince us 
that the emission band B is due to a phonon-assisted process, we calculated the sum of occupation numbers 
for both equilibrium phonons (lattice temperature T[J and nonequilibrium phonons (carrier temperature T c ) by 
assuming the carrier temperature is the same as the effective temperature for the nonequilibrium phonons 
after 30 p3 (see later discussion on the carrier temperature). The calculated result is shown by the solid 
line in Fig. 3. The fitting of the total phonon occupation number to the obtained intensity ratio from the 
spectra is impressively good. This implies that the intensity of peak B is well correlated with the L0- 
phonon population. Because of participation of nonequilibrium phonons the intensity of the B band (see also 
time-resolved luminescence at 780 nm in Fig. 4) is strongly enhanced especially at low temperatures. This 
further supports our assignment of the B band to the nonequilibrium phonon enhanced phonon replica. 


(ill) It is not possible to attribute the emission band B to the emission from the renormalized band-band 
transition because the peak positions of the A band (766.7nm) and B band as well as the energy separation 
between them at 4.3K do not change with variation of the excitation intensity by a factor of 9x10”' (2ND3). 
Further support to the above statement is the fact that the spectral position of peak-A exactly coincides 
with each other for two different luminescence studies using different light excitation sources; the one is 
the weak 200 ts-pulse train with a repetition rate of 114 MHz, while the other is the amplified 0.5 ps-pulse 
with a 20Hz repetition rate. As the observed band edge transition energy does not change over a range of 
photogenerated carrier densities (10 l *-10 , Vcm*) the bandgap renormalization in this study must be excluded. 

(Iv) We have also measured the polarization of the emission bands A, B, and C by adopting a "right- 
angle" ls >“ geometry for obtaining integrated luminescence spectra. It was found that the emission of band B 
detected along the y-direction (emitted at the sample edge) was strongly polarized in the x-direction, while 
the A and C emission bands are depolarized independent of excitation power density and lattice temperature. 
The intensity ratio between the x- and z-polarization for the B band was about 20 and independent of lattice 
temperatures from 4.3 to 300K at full excitation power density P m . But it is strongly dependent on the 
excitation power density at a given lattice temperature. For example at T L -100K, the B band cannot be 
clearly identified from the A band and shows no difference between the two polarizations at the excitation 
power density 0.017 P m (2ND13), while the ratio at full excitation P m is 20. The polarization behavior of the 
B band is consistent with the Raman scattering measurements reported by Zucker et al." where the initial 
photons are provided by the external laser source. In our case, the initial photons are provided by the 
luminescence from the recombination of n-1 electrons and n-1 heavy-holes at the subband edges. This process 
is 3trongly enhanced by the presence of a large number of nonequilibrium phonons emitted by the hot- 
electrons at high excitation. The polarization behavior of emission band C is expected since the emission is 
from the 1.2pm GaAs buffer layer. This, in fact, provides a good check of our system. However, the 
polarization behavior of the emission band A is not understood in terms of the selection rules of dipole 
recombination described by Iwamnra et al. 1T , where the emission intensity for x-polarization (TE polarization) 
should be at least four times larger than that for z-polarization (TM polarization) due to different density 
of states for heavy-holes and light-hole3. 
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There has been much debate 1 ’" 1 ’ on the interpretation of the spectral feature below the n-1 electron 
heavy-hole transition energy. This is because the energy levels for Impurity states, longitudinal optical 
phonon (LO)-phonon replica, and intrinsic bandgap-reduction due to many-body interactions roughly coincide 
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uith each other. Holonyak and co-workers'* have demonstrated the LO-phonon participation In GaAs-AlGaAs 
based QW laser emission by the observation of more than one LO-phonon sideband of laser operation below the 
n*1 confined-particle transition. Also, the phonon replica of n*l electron to heavy-hole emission has been 
observed in the InGaAs based systems. The observation of the coupling at both the GaAs- and Ir.As-like bulk 
LO-phonon energies of InGaAs QW by Skolnick et al. 20 has confirmed that the peak below the n = l electron- 
heavy-hole transition is a phonon replica. Furthermore, it should be pointed out that our assignment for the 
B emission band is consistent with the time-resolved pnotoluminescence data described in the next section. 

b) Time-Resolved Luminescence 

Time-resolved photoluminescence emitted from the symmetric MQW structure at 3K are shown in Fig. U for 
various emitted photon energies. The emission centered at 770 nm arises from the recombination of electrons 
and heavy-holes near the respective first subband-edges. The emission with wavelengths less than 770 nm is 
from the recombinations of hot carriers. The time-resolved luminescence intensity assigned to the phonon 
replica 3t 780 nm is shown in Fig. k. Each luminescence profile was the average of 10 individual shots using 
a prepulse 2 ' for averaging. The left peak of the dotted curve is the prepulse which reflects the 2 ps 
temporal resolution. The right peak on the same curve is the Rayleigh scattering light from the sample 
surface which defines the "zero" time for our analysis. Several features appear in the data displayed in Fig. 

(i) Although the risetime for all the luminescence profiles is instrumental, the 2 ps uppei—limit 
reflects the rapid thermalization and initial cooling. (ii) The shape of the rise-time of the luminescence 
profile at 780 nm is similar to that at 770 nm, but delayed by 3 ps. This implies that the emission at 780 
nm does not originate from the same band as the emission at 770 nm; otherwise, the rise-time of the 
luminescence at 780 should start at the same zero point as that at 770 nm. The 3 ps delay is consistent 
with our assignment of the emission band center at 790 nm to the phonon band. This time is just the 
effective time required to establish a large nonequilibrium hot phonon population. (iii) The luminescence 
decay time of 30 ps (broken curve) for the phonon replica emission at 780 nm cannot be due to impurity 
emission at high excitation. 

The decay times were obtained by fitting the data in Fig. a to the expression given by: 

I(E,t)-I,[l-e- 1/Tr je t/Td , ( 1 ) 
where the parameters (I 0 ,t r T d^ are the proportional constant, the risetime, and the decay time, respectively. 

To study the energy relaxation process quantitatively, the earner temperature as a function of time must 
be determined. There are many ways to obtain the time evolution of carrier temperature. One way is to 
analyze the high energy tail of time-resolved spectra using a Maxwel 1-Boltzmann distribution. 6- * Another way 
is to measure the ratio of luminescence intensity profiles as a function or time at two energy positions 
within the exponential energy tail by also assuming a Maxwell-Boltzmann distribution. The above two methods 
cannot extract information on the photogenerated carrier density. An expression in the time domain is 
introduced to fit the experimental luminescence temporal profiles by using two adjustable parameters, 
namely, the carrier densities (n e , npp, n^p) and the electron temperature (T c ). 

Assuming the direct optical transition, the luminescence intensity is given by: 

I( Ei,t) =Cj( 1-e r )[|M e -pp i*Pe ^e Phh ^ h * l^e-lhl J °e ?e Plh ^h^* (2) 

where p e pp lh " m e hh lh is the density of state for electrons, heavy-holes, and light-holes; M p -nh(lh) is 
the matrix 'element 'for 2 electron to hh(lh) transitions; Cj absorbs all the constant factors including the 
corrections for detector response and the transmission of each narrow band filter used; t r is the risetime of 
luminescence and is set to be 1 ps (approximately the response time of the system) which is an upper-limit 
of thermalization time of the electron-hole system; 

r e p - i/( e te i" u e,h> / K T c . i) (3) 

Is the Ferml-Dirac distribution for electrons in the conduction band (with subscript e) and for holes in the 
valence band (with subscript h). The quasi-Fermi energy of holes (up) is related to the quasi-Fermi energy 
of electrons (u^) and the carrier temperature by the relation. 


n e' n hh * n lh 


This value 13 given by: 


u h - KT c ln[e aln <’*e Ue/KT °)-ll 


An unique set of parameters of T c (t) and u e (t) was used to consistently fit all the luminescence profiles 
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detected at different photon energies Ej. Three typical calculated luminescence profiles are shown in Fig. 14 
by the thick solid curves. 

In Fig. 5 the time-resolved photoluminescence profiles from the asymmetric quantum well are shown at 
three typical energies. The emission centered at 760 nn arises from the recombinations of electrons and 
heavy-holes near the respective first subband-edges. The theoretical solid curves in Fig. 5 were generated 
using Eq.C2). The same set of values for T c and u p used to fit the data for the symmetric quantum well were 
used to fit the luminescence profiles of the asymmetric quantum well. The energy relaxation process is 
similar for both wells despite the different shape of the well. This is as expected! 

The experimentally determined values of T c as a function of time are plotted as the solid curve in Fig. 6 . 
The shaded area reflects the extent of uncertainty in deducing the cannier temperature within the first H ps 
due to our limited time resolution. The plot of T c vs t is interpreted as follows: carrier-carrier 

collisions quickly lead to a th®rnaliz®d distribution at very nigh T c within 1 ps due to a large number of 
exeitec carriers . 22 The thermalization process of carriers can only be probed by femtosecond 

spectroscopy. The Initial cooling of the thermalized distribution is studied with our present time 
resolution of 2 ps, providing information about both the electron-phonon interaction and the carrier 
lifetime. The initial cooling within tne first ‘3 ps is very fast (25CK/ps) assuming a Fermi-Dirac 
distribution for electrons. Aft®r the initial cooling, a large number of longitudinal optical phonons in 
finite wavevector spac® ar® ucoumul it®d u- to the finite pr.onon lifetime . 21 The phonons emitted by the hot 
electrons may be reabsorbed by the electrons as the reverse process of emission giving rise to a slower 
cooling for the hot-earn 1-rs. It sr.oul 5 :•> -mpr. isized r.®r® that the rapid initial cooling rules out the 


importance of significant scr®®r..-.g of tne -l-ctroi 
screening was important, the initial -an:: ceding 
not te built up. Another important Observation Ais 
slow component of the carrier temper >* jr- cooling 
pr.onon enhanced phonon replica emission at 730 cm 
and tne nonequilibnium phor.or, systems to -quilisr. 


decay constant of 30 


tne el-ctron-pr.onon interaction in the present study. If the initial 
inc -xdsng would be slow and a nonequilibrium pr.onon population would 
servati :r. displayed ir. Fig. 0 is that the time constant - 30 ps for the 
' ur- cooling curve is the same as the decay time of the nonequilibrium 
:n at 700 nm. This suggests that tne time required for the electron 
to equilibrate with earn other is k-3 PS. The system then decays with 
coupled hot ®lectron-tnonon system. 


Th® experimentally determined quasi-Fermi-er.ergies for electrons and holes are plotted in Fig. 7. The 
changes of u p and u^ are very rapid within th® first 10 ps and the behavior of the degeneracy for electrons 
and holes as a function of time ar® reversed. 


Using the experimentally determined distritutions f 
the energy loss rates for electrons and holes as functi 
is defined as follows: 


electrons and holes in terms of T c (t) and u Pifl (t) 
>f T c or a,, can be obtained. The energy loss rate 


f (T c 

• u e,hh, c 


However, if the electron-hole system is treated as a Maxwell-Boltzmann gas, i.e., the electrons and holes 
have the same distribution function which is a Maxwell-Boltzmann distribution, the energy loss rate for the 
carriers (c-e,hh) is given by: 


It should be pointed out that the energy loss rates for electrons and heavy-holes obtained by Eq. (7) are 
the net loss rate3. It includes the energy loss due to electron(hoi®)-phonon interactions and the energy 
exchange between electrons and holes by carriei—carrier scattering. The hole-phonon scattering rates were 
calculated 2 ' and measured 1 tc be 2.5 to 3 times larger than the rate for electrons due to the additional 
coupling through the deformation-potential. Therefore, an energy transfer process from electrons to holes 
should be expected in our photogenerated carrier system. The experimentally determined energy loss rates 
using Eq. (7) or Eq. (8) are shown in Fig. 8. The solid and dotted curves were obtained from Eq. (7) for 
electrons and heavy-holes, respectively. The dot-dashed curve was obtained from Eq. ( 8 ). It should be 
pointed out that the curves obtained from either Eq. (7) or Eq. (8) are the experimental data because the 
distribution f(T c ,u Pi h) was experimentally determined. The broken line is calculated based on the simple 
theory 1 * where lattice can be treated a3 a heat bath for quasi-equilibrium carriers. This theoretical curve 
does not agree with experimentally determined curves. The smaller P hh and P p is because of the presence of 
light-hole population and the different masses for electron and hh which lead to different behaviors of 
qua3i-Fermi energies for electrons and heavy hole3 as a function of time. 
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The salient features of the experimental data displayed in Fig. 8 are: 

(i) The experimentally obtained energy loss rate for electrons is more than two orders smaller than 
predicted by the simple theory 5 (broken line). The reasons for thp differences arise from the reabsorption 
of nonequilibrium phonons by electrons. As stated above, the energy loss rate determined for electrons by 
scattering alone with LO-phonons via the wave-vector dependent Frohlich interaction may be even smaller than 
our determination because the net energy loss rate determined from Eq. (7) is the sum of rates caused by 
electron-phonon and electron-hole interactions. 

(ii) The net energy loss rates for electrons and heavy-holes almost follow the rate obtained from Eq. (8). 
The reasons for this are the rapid decrease of u p as shown in Fig. (8) and the slow cooling of the carrier 
temperature which results in u e /KT c < 0.5. This makes the electron-hole system behave more like a Boltzmann 
gas. 

(iii) The energy loss rates increase as T c increase. The dot-dashed curve even bends over and approaches 
the maximum value of 5.3xlO"‘W for bulk GaAs. This "bending” behavior reflects the initial rapid cooling of 
hot carriers and was also observed by Shah et al. 2 This may be the indication for the presence of 
nonequilibrium phonons in the system. 


In Fig. 9, we have plotted the average phonon emission time for hot carriers defined by:’ 

T . Em e - E Lo /KT c 

T bVg - p 


(9) 


as a function of the carrier temperature. This is just another way to describe the energy relaxation 
process. As can be seen the ^avg i 3 carrier temperature dependent and thprpfore time dependent. For 
T C >1 200K or t<5 ps the T av e remains at a value of about ! ps. As time t>5 ps (or T c <12CCK) the T av ~ 
increases quickly with decrease of the carrier temperature. At T c = 980K, i 2V g is 10 P 3 - Ryan et al. 
obtained a constant value of 7 ps for t av? by studies of time-resolved photoluminescence of modulation 
doped-MCW. The reason for the difference is Simply because their time resolution did not allow them to 
monitor the entire carrier cooling process. Our upper limit of 0.9 ps for the initial carrier cooling agrees 
with the theoretically predicted’ value of 0.16 ps within our resolution. As soon as a large 

nonequilibrium phonon population is built up the carrier cooling is suppressed. 


We have extensively discussed the energy relaxation process based on the experimentally determined 
distribution function f e p(T c ,u e ,e ). In the same manner information about the carrier lifetime can be also 
obtained using determined f e ft. ' In Fig. 10 th? carrier densities for electrons, heavy-holes, and light-roles 
as a function of time are' determined from the experimental data. The densities are obtained by the 
following expressions: 


n e (t) = p e f e d c (10) 

. 0 

"hh^ * Phh f h d e (11) 

Jo 

n !h (t) * p lh r h d e (12) 

. AE 

where AE is the energy separation between the hh- and lh-subbands at the zone center. The relation n e * npp 
♦ n^h must be satisfied. The quantity of In [n e (t)j is also plotted in Fig. 13 as a broken curve. 

The salient feature of the data in Fig. 10 is that the carrier density decreases nonexponentially and very 
rapidly within the first 30 ps after the end of a 0.5 ps pulse excitation. One cannot define a single carrier 
lifetime due to the nonexponential nature of the density decay curves. However, an effective carrier 
depletion time (time for the density to decrease by a factor of e _l from n e (t-0)) can be deduced to be as 

3hort as 10 ps. This decrease of carrier density occurring in such a short time cannot be accounted for by 

the usual blmolecular recombination, which is a much slower process, on the nanosecond time scale in the 
bulk GaAs and about 350 ps In quantum well structure with well thickness of - 5 nm.” This short carrier 
depletion time Is not caused by the rapid carrier diffusion which will be discussed in the diffusion section. 

Such a short carrier depletion time can be explained consistently by the participation of nonequilibrium 
phonons. The nonequilibrium L0 phonons In the quantum well system can play two essential roles. First, 
they substantially suppress the en rgy relaxation process of carriers due to reabsorption of these phonons 
by carriers. Second, since these phonons are localized in both wavevt"tor space and the real physical space 

(In the well plane) they will behave like coherent Bosons. The existence of the coherent Bosons will tend to 
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increase the number of the bosons in the system with a rate proportional to the present number of Bosons 
causing a rapid carrier relaxation. This is a stimulated phonon emission process and can only be 
accomplished by the phonon replica emission, i.e., when an electron in the conduction subband recombines with 
a hole in the valence subband not only a photon but also a phonon or more phonons will be emitted to Join 
the Boson system. Since this is a stimulated process, it is very effective in exhausting carriers from the 
system. Unlike the conventional stimulated emission in a semiconductor laser, this process is expected to be 
even more effective at room temperature since equilibrium phonon occupation will also increase as lattice 
temperature increases. We have confirmed the nonequilibrium phonon induced stimulated-phonon-emission by the 
detailed polarization studies of integrated luminescence spectra discussed in the last section. By employing 
the nonequilibrium-phonon model we are able to explain not only the slow carrier cooling, but also the 
ultrashort carrier depletion time in highly photoexcited undoped quantum well structures. 

Carrier Diffusion 

Electron-hole plasma expansion at high photoexcitation power density in semiconductor bulk and quantum 
well structures is a subject of great interest from both the points of view of theory and device application. 
Various methods have been employed to study the dynamic process of carrier diffusion such as spatially and 
spectrally time-resolved luminescence spectroscopy* 4 - 27 and Raman scattering. 11 - 2 * Different mechanisms have 
been suggested to explain the rapid carrier diffusion in bulk semiconductors. Rapid electron-hole diffusion 
has been observed 27 in CdSe at RK but not at room temperature. The absence of a larger number of thermal 
phonons at low temperature made it possible to observe the rapid carrier diffusion process. Using a streak 
camera, Junnarkar and Alfano 2 ’ were able to measure the difference of spatial profiles of the pump laser and 
the photoluminescence at various excitation levels by which the diffusion velocity was estimated. Junnarkar 
and Alfano 2 ’ obtained a diffusion velocity on the order of 10*cm/s at high photoexcitation even at room 
temperature. Such a large diffusion velocity was attributed to screening of polar electron-phonon 
interactions and partially screened nonpolar electron-phonon interactions. High internal pressure of a very 
dense electron-hole plasma has been proposed 20 to account for the rapid expansion of plasma. The pioneering 
work of Tsen and Morkoc 11 using the time-resolved Raman scattering technique has shown no Influence on 
photogenerated carrier diffusion by carrier confinement, zone folding of phonon modes, and phonon confinement 
in symmetric CaAs/AlGaAs MQW structure at low carrier density. In view of their works we should point out 
that the well thickness of their MQW structure was so thick that any confinement effects would be smeared 
out; and the carrier density deduced from Raman data was the total carrier density within whole energy 
distribution. Therefore, they could not obtain the information about dynamic diffusion of carriers with 
different kinetic energies. In the following subsection the results on dynamics of carrier diffusion 
processes in the x-y plane (lateral) will be presented. 

A schematic diagram of the photoluminescence image on the screen of the streak camera is depicted in Fig. 
11 to show how we study dynamics of carrier diffusion. The area inside the filled-llne circle is the laser 
excited area. If there is no carrier diffusion photoluminescence should be restricted to inside this circle. 
Suppose sufficient electrons and holes can diffuse far enougn to the unexcited area where they can recombine 
one should be able to detect the luminescence from the unexcited area (outside the filled-line circle). 
Using the streak camera we were able to measure the laser spatial profile and luminescence profile at the 
sample surface. Those spat.ul intensity profiles were integrated over a time duration of 1 ns. The width 
(FWHM) of spatial profile should give information on the carrier diffusion length within 1 ns. However, this 
information is not adequate to deduce accurate values of diffusion constant and diffusion velocity. To study 
the dynamic diffusion process quantitatively we make use of two electronic windows to obtain time-resolved 
photoluminescence from the excited (I 3C ) and unexcited (I se ) areas at the same time (and thus laser intensity 
fluctuation was completely avoided). 

The time-integrated laser and luminescence profiles at the surface of the 5bi symmetric MQW structure are 
plotted in Fig. 12. The dotted-curve is the laser profile. The inside and outside solid curves are the 
luminescence profiles at 7A0 nm and 770 nm, respectively. The broken-curve is the luminescence profile at 
830 nm, which comes from the 1.2 urn buffer layer. The gain of the detecting system was adjusted such that 
the maximum Intensity recorded was roughly the same. The widths of luminescence profiles above the subband 
edge (760-700 nm) are exactly the same as the width of the laser profile (dotted-line) reflecting no or 
little diffusion for high energy carriers within our spatial resolution. The width of the luminescence 
profile at the subband edge (770 nm) is about 2 times wider than the width of the laser profile. This energy 
dependence of the carrier diffusion process is not usually expected because a larger kinetic energy is 
associated with larger carrier velocities. However, the nonequilibrium phonons in the system during the 
first 30 ps after the laser excitation may inhibit carrier diffusion, especially for the carriers with high 
energy. The extent of carrier diffusion at the subband edge seems to be consistent with the results given by 
Tsen and Morkoc. 11 The diffusion velocity determined by them 11 is 8.6x10‘cm/s. According to this value the 
carriers will travel about 100 pm within 1 ns which will broaden the width or luminescence profile from 200 
tim to about A00 urn. This agrees with the measured lateral width or the luminescence profile at 770 nm in 
Fig. 12. It should be pointed out that the Initial carrier density in our experiment was about 100 times 
larger than in their experiment. Due to the short carrier depletion time (10 ps) in our case, the average 
carrier density within 1 ns would be similar in both cases. 

The time-integrated laser and luminescence profiles at the surface or the asymmetric quantum well are 
plotted in Fig. 1 3. The Innermost solid curve and dotted curve are the laser profile and the Caussian-line- 
flt of the laser profile, respectively. The center wavelengths of the luminescence are indicated on the 
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curves. The key differences of the spatial luminescence intensity profiles of the asymmetric well from that 
of the symmetric wells are: the widths of luminescence profiles are much wider than the width of the laser 
profile; and the width of the luminescence profiles do not depend much on the wavelength of luminescence. 
The possiole mechanisms are given below. 

In Fig. Id, we have plotted the time-resolved photoluminescence from I se and I sc measured at 4K with 2 ps 
time resolution for asymmetric well. Based on the fact that the width of the spatial luminescence profile 
does not depend on the photon energy detected within our spatial resolution, w» placed a wide-band filter 
centered at 700 nm in the luminescence path to increase the strength of the luminescence intensity. The 
broken curve is the Rayleigh scattering light from the sample surface. As can be seen, the risetime of the 
I 3C is time-resolution limited, while the rise of I se is much slower than that of I s= 3 nd also well resolved. 
The fast decay of I so has been discussed in the energy relaxation section. It should Be emphasized here that 
using the same set of experimentally determined T c (t) and u e (t) we have determined the carrier depletion time 
is 10 ps. The only choice for theoretical fittings to I se is the use of stancard classical diffusion equation 
to calculate the carrier density as a function of time t and distance r. The exact values of separation 
between the two windows and the widths of the windows are given in Fig. ' I. This method should be adequate 
as long as electrons act more like particles rather than waves. 

The diffusion equation is given by: 


where t is the carrier lifetime, and D is the diffusion constant of carriers, 
carrier density is assumed 27 to be the same as the laser profile which is: 


The initial condition for 


n(r,t= 0 ) - n 0 e J 1 

where n 0 is the carrier density at the center of the laser spot at t=0, ind a is 1 ZZ pm. ct 
the laser profile to Eq. (I 1 !). The solution of the second order differential Eq. {'■'}) with 
condition Eq. (14) is given by : 27 


n(r,t) - n 0 e " 7 l 


air.ed by fitting 
he given initial 






^ . I 6 lcg 7 Dtj-i/ 

1 r‘/( a 2 / J log 2 * a Dt) 

(15) 

i 

a 


Assuming ambipolar diffusion, i.e. electrons and holes have tne same diffus an velocities, 


Using Eq. (1 6 ) the parameters of D and t have teen determined to be 10*cm 2 /s and 11 ps, respectively. 

The results are very surprising in two respects. First, the value fcr 0 is more than four orders of 
magnitude greater than the conventional diffusivity. 71 ° diffusion velocity v ■ »\r/ T ) . i.JxlO’cm/s is four 
times larger than the initial Fermi-velo-ity. It shoulJ Be pointed out that the values for D and t deduced 
from the time-resolved luminescence profile I se with the theoretical expression Eq. (lb) are consistent with 
time-integrated spatial luminescence profile. Second, the carrier lifetime at t<25 ps deduced from the 
diffusion Eq. (13) i 3 the same as the carrier depletion time determined from the time-resolved luminescence 
data shown in Fig. 4. Thl3 demonstrates that our analysis is sel r-consistent. Using the determined 
n(r,t) we can show the contribution of the second term m the right-side of diffusion Eq. ( 13 ) in losing the 
carrier at a given position r is always smaller than the first tern and thus the ultrashort carrier depletion 
time cannot be attributed to the carrier diffusion tut to the nonequilibrium phonon ennanced phonon replica 
emission. 


Several possible origins of the extraordinary rapid carrier diffusion in the x -y plane of the asymmetric 
quantum well can be given at this stage. (1) The internal pressure 10 of high electron-hole pairs may be 
extremely large when the density is well 3 bove the electron-hole liquid density. (2) There are a large 
number of the coherent nonequilibrium phonons emitted by hot electrons in the system. These phonons have 
momenta parallel to the x-y plane and may also drive the carriers outside the excited region by adding the 
phonon momenta to the carriers. (3) Due to the asymmetric potential well profile there are no transmission 
states (virtual states or unconfined states) above the energy barriers unlike the symmetric well case. This 
la the intrinsic property of asymmetric quantum wells. The photogenerated carriers may be Immediately 
scattered by the higher energy barrier and obtain the momentum from the interface which assists carriers to 
expand in the x-y plane. Due to a large kinetic energy or -300 meV, the photogenerated carriers acquire a 
la ge initial diffusion velocity. The origins (1) and (2) seem to be inconsistent with the symmetric well 
case. However, whether the mechanism described in (3) is responsible Tor the observed rapid diffusion is 
still an open question worthy of further study. 
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Conclusion 


In this paper extensive data have been presented on time-resolved photoluminescence studies as well as 
other related investigations of a symmetric MQW structure and an asymmetric quantum well with similar well 
thickness. With 2 ps time resolution we were able to study the initial energy relaxation and carrier 
diffusion processes simultaneously under the same conditions. 

The existence of a large population of nonequilibrium phonons in highly excited semiconductor quantum well 
structures is experimentally verified. Its effect on the energy relaxation is to slew down the cooling rate 
after an initial rapid cooling (0-5 ps). The symmetric and asymmetric potential well profiles do not show 
any difference in the energy relaxation process under the same conditions. 

A new mechanism to explain the ultrashort carrier depletion time deduced from the consistent fitting of 
time-resolved photoluminescence profiles at various emitted photon energies, is proposed to be nonequilibrium 
phonon enhanced phonon replica emission rather tnan other nonradiative processes suen as Auger recombination, 
carrier diffusion, and the conventional radiative stimulated emission. 

The differences in carrier diffusion for two different potential well shapes are remarkable. For the 
first time we have observed the ultrafast carrier diffusion in the x-y plane of the asymmetric quantum well, 
while there is little or no carrier diffusion except for carriers at the subband edge in the x-y plane of the 
symmetric quantum well under similar experimental conditions. 
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non-commutation relationship between the Hamiltonian of the system and the 
operator in the confinement --direction was taken into account using random variables In the calculation of 

electron scattering by polar optical phonons in quantum wells. The rate at the onset of phonon emission is 

found to be 3 times smaller than previously reported. Furthermore, this rate becomes on increasing function 

of initial electron energy instead of a decreasing function. 

It is important to take into account the effect of electron confinement on the scattering of electrons by 
polar optical phonons (PC?) in quasi-two-dimensional (Q2D) quantum wells. 1 '” The Hamiltonian H of the 
electron system and the oper itor k~ for the electron momentum in z-dmection do not commute with each ''then 
because of the potential profile experienced by electrons in the confinement z-directior.. That is, 


electron momentum 
Che calculation of 
phonon emission is 
ncreasing function 


he usual conservation law 


momentum in z-direct ion: 


must be modified m the sense that q z , k'- and k, are r.o longer deterministic observable variables but are 
random var.abies .r.v., 1 . here, q^ is the z-ccnponen' of tne pr.on-n wave vector, and k'-. are 'he electron 

wave vectors tefore and after scattering, respectively. The statistics of each r.v. is characterized by its 
probability density function. The pioneering wore of Hess 1 , Price 2 and Leourton” essentially adopted a "non- 
m omentum-select ion rule" in z-direction for calculation of phonon emission and phonon absorption rates. They 
-reat c^ as 5 'ip de t - rr.i.i j v ?. ►;*. A1 1 ,jrr. ise i y J poaiio ; :i " fj;• y t or. ?, j— censervst ion" in z~ 

direction, he used instead a "momentum conservation approximation". Furthermore, the value of k, in wave 
function in the coordinate representation was eor.zi ~.vr>-o is tne -iectror. momentum m z-direction in his 
analysis. This approach is r.ot valid since < z is equal to a prv-ise value of n/L, in an "infinite well" 
approximation and for a "finite well" g z is just a parameter tr.at related to the eigen-value of the 
Hamiltonian of tne electron system. 


Furtr.ormcro, the value 


by its 


a "non- 


They 

«L 

n” :n z- 


in wave 


m his 

e well" 

/'v““ 


The thesis of this paper is to present a theoretical analysis of electron scattering by PCP in Q2D 
structures placing emphasis on confinement effects by ta-'ing into account the non-commutation relationship 
between the Hamiltonian -per s' or H ar.d the electron momentum opera*-p k.. We show now "fuzzy mementum 
conservation" in z-direction is exactly feateq using rcitner the "nor.-momer.tum-seiection rule" nor the 
"momentum conservation approximation" but by consider.ng as i random variable. A. key result that 

distir.guisnes this work from tne previous works 1- ” is that the rate at the onset of pnonon. emission is 
reduced by a factor of 3. 


There are two important aspects of tne non-s'mmut a-ion reiitmshir given by Eq. (1). First, the 
corresponding observables cannot be measured simultaneously with arbitrary accuracy. The product of their 
quantum root-mean-squ.are deviation must rot be less than tne naif of tr.e quantum expectation value of the 
dV(z) 

operator i in a given state. Second, the electron energy Vue to the motion in z-direction has no 

correlation with the corresponding electron momenta- in this direction. The eigen-value of momentum 

operator k z cannot precisely be determined due to the electron confinement by the potential well and must be 
considered as a random- variable (r.v.). When an electron is restricted in 'he first bound state the quantum 

expectation value ;f lK,w zJ is zero resulting in „ precise determination of the lowest eigenvalue of operator 

H and root-mean-square deviation for k z car, te deter med. Those two points play the key role in 
calculating the electrcn-POP scattering rate. 

In our analysis, the unscreened electron-PCP interaction, equilibrium bulk phonon spectrum and an infinite 
depth square well with the lowest subband occupied are assumed. 

Within the framework of the effeetlve-mass- jpproion, the wave function for electron in coordinate 
representation is then given by 

i pi kr 

|k d(z). (,) 

/A 
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is the wave function in z-direction describing the bound state. The vectors r=(x,y) and k = (k x ,ky) denote the 
two dimensional electron position and the wavevector parallel to the interfaces, respectively. The 
corresponding electron energy can be expressed as 


h k 
2m 


(5) 


E z-irT-,. (6) 

z 2mL z 2 

Here m is the isotropic electron effective mass, L z is the well width and A is the layer area. 

The probability density function for the r.v. k z can be obtained from the wave function given by Eq. (U) by 
a representation change to the wavevector space using 


(7) 


$(k 7 ) =- o(z)e" 

/2 77 1 -o= 


This yields 


p(k z ) = |<|>(k z )| 


os 2 (k-L, 


ii (1 — k z L z 


77 2 ) 2 


( 8 ) 


Since there exists a "fuzzy momentum conservation law"' in z-direction given by Eq. (2), the probability 
density function for the r.v. q z , p(q z ), can be obtained from the convolution expression: 


p(q z ) * p( ik z ) * p(*k’ 2 ). 


( 9 ) 


The electron scattering rate (w k . k ) per unit area in k-space in 52D structures can be expressed in Born 
approximation by Fermi's golden rule by: 


h(2 


TP I M x’x | 26(El < " 


- h ui q) > 


'C(nn*1)‘ 




UJ 


( 10 ) 


( 11 ) 


( 12 ) 


is the matrix element for the Frohlich electron-POP interaction. For simplicity we have only considered 
phonon emission process. Q and V s denote the three dimensional phonon wavevector (q,q z ) and the sample 
volume, respectively. The 6 is a modified Kronecker function defined as 


ll if a*b 
(0 if a-b 


f.notion F in Eq. (11) is the form factor in z-direction which is defined as: 


F^dz’^z) “ ~ 
L z 


-iq z z 


sin -— dz. 
L z 


(13) 


(19) 


re is to realize that the matrix element M k > k is the function of a r.v. q z . 
c , r i v-j f 793 U/tra/ast Lasers Probe Phenomena in Bulk and Mtcrostruclure Semiconductors (19871 / 71 


\\V, 
'«• "w 1 , 
f.-h. 


wJ&Sj 

M 

4 


1 

m 


' • 

r>:% 

v\S 

'-''M 

AA 


-> s 


m 

yy 

► V V 

cSa 


ll 













Therefore, the w k > k also becomes a r.v.. The most likely value of the w k . k that can be experimentally 
measured is given by its expected value using the statistics associated with the r.v. q,. This is. 


w k'k * I P(Qz) w k'k 


The total emission rate is then given by 


"k'k = i w k’k d q , I^~\q, 1^ y £{ E k -E k * h ujo)z k i_ k ,q qdq 


lQ 2D (q,L z ) = I 0 


and 0 denotes the angle between the k and q. The factor I„ is defined 3S: 


The function g(x,y) is given by: 


, sin 2 (—) 

g<x,y ■ 0 - 1 75775 (f -7- 09 ) 

9 x j 1"(xz / i«) , ! (xz / 2) 2 

cos 2 j~(x*y)z /2~jcos 2 (yz /2) 

- (x *y ) 2 z 2 /^j^l-y^ 2 / n 2 J 2 

with Z“q„L z , which includes 1 / Q 2 from Frohlich electrcn-POP interaction Hamiltonian, the form factor in z- 
direction and the probability density function p(q z ) f or r.v. q z . All the wave vector are normalized by 
/2m uiq 

q 0 *- - L„ - ‘ which is about 1 / AoU for GaAs. 5 

h z 

The integral lQ 2D given by Eq. (17) for Q2D structures represents the electron-PGP interaction strength as 
function of q. It contains all the information about confinement effects as well as size effects' on the 
scattering rates. In order to further understand how electron confinement makes the difference in the 
electron scattering rate between Q2D structures and its bulk counterpart, a plot of the electron-POP 

interaction strength l3D versus q for 3D electron system is needed in such a form that can be compared with 
I 2D . It can be shown that 

I„ I„ 

j3D(q) , —-— * — ——, (20) 

Q + 2 . ^ ^ Q * Ak 2 ( ) 

where Ak z (')-(k z *k’ z ). The first and second terms on the right side of Eq. (20) represent the scattering 
processes without and with change of initial direction of electron momentum, respectively. Since Ak,(*) in 
Eq. (20) depends on bojh the magnitudes as well as the scattering angle 6, one cannot find an explicit 

expression for for Ak z ^ ) in terms of q and 8 unless a variational parameter for the electron kinetic energy 
ratio (r- K 2 / k z 2 ) between the total kinetic energy and the energy due to the z-direction motion is 
introduced. 

Using the momentum and the energy conservation laws for the 3D scattering process, the Ak z ( + ) can be 

expressed as 

- (q 2 *0* /(q 2 *C) 2 -f,Mq 2 cos 2 0 

ak z ( 1 " teet- (2D 

2qcos 0 / r-1 


where OSOSarccos 


Vir and 


£" r ~ r w ith 1 < r < ®. The value of r = 3 implies that the kinetic electron 
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energy is equally distributed in all directions, while r<3 more electron kinetic energy is in the z-direction. 
The 2D electron motion in x-y plane corresponds to the case when r approaches infinite. 


We now present the numerical results in the following sections. 


The two solid curves in Fig. 1 show the probability density profiles for the r.v. q z obtained from Eq. (9) 
for well width L z * 2L 0 and 5L 0 . The larger well width L z , the smaller is the variance of the density curve 
resulting in the less fuzzy momentum conservation in z-direction. The magnitudes of F(q z ,L z ) versus q, for 
L z = 2L 0 and 5L 0 are also displayed in Fig. 1 as two dotted curves. The two features from the plots should 
be noticed. First, in the limit q z -*0, |F(q z ,L z )| approaches unity independent of L z reflecting a scattering 
event with an exact momentum conservation in z-direction. This is also the most favorable process according 
to the density curve for q z . Second, for very thick well the |F(q z ,L z )| can be written as 6 k „_ k i, <, which 
implies the momentum conservation law for any scattering event and p(q z ) approaches a 6-furfctioTa ?or the 
bulk situation. The fact that not every scattering event follows the exact momentum conservation law in z- 
directron is due to the electron confinement. This procuces the weakening of electron-PCP interaction in Q2D 
structures. 


The solid curves in Fig. 2 show the integral lQ 2D as function of the parallel component of phonon wave 
vector q with L z as a parameter. Curves 1, 2, 3, and k are plotted in Fig. 2 for the value of L- of C.5L„, 
L 0 , 3L„ and 5L 0 , respectively. Our calculation procedure using the probability density function p(qO for r.v. 
q z yields a significant reduction by a factor 2 on the value of the integral I'w-D obtained by Leburtcn" in the 
range from 0.1 to k for q. In spite of this key difference other features of our integral T' 2D are very 
similar to what Leburton has shown." Namely, the integral lG 2D is a monotonic decreasing function of q for 
Q2D electron system and weakly dependent on L z . This similarity occurs just because the profile of 
probability density function for q z happens to have similar shape as |F(q,,L z 5|. 


The interaction strengths l3D f or S y Stem with various values of 8 and r are also displayed in Fig. 2 by 
the dotted curves 5 to 8. The values of r for curves 5, 6, and 8 are 10, 3, and 1.2, respectively. The 
value of 0 for those three curves is set to be zero. The curve 7 with r = '.2 and 6 = it/ 6 is used to see the 
variation of l3 D (q) by the change of the scattering angle 8. A comparison of the dotted curves 5 to 3 for 3D 
with the solid curves for Q2D provides the physical insight into electron confinement in quantum wells. The 
following three points are discussea: first, the value of 1-3° increases with increase of the value of r. It 
indicates that the electron-POP interaction strength for given value of q in the system with r>3 is stronger 
than the system with r<3. The reason for this is the-weawness of electron motion m z-direction resulting in 
small amount of change of the electron momentum Ak z ( K It should be emphasized that the value of T^D a t q 
as L z ->0 do not approach the value of l3D with r-*« which corresponds to tne _D case. In other words, as far 
as the electron scattering process is concerned the C2D quantum well structure with L z -*0 is not equivalent 
to a 2D system due to the potential barriers which confine electrons inside the well. The values of and 
I<32D are close to each other for the value of the ratio r between 1.2 to l.’i. This implies the origin of the 
difference in electron scattering between Q2D and 3D systems comes from tne fact that the electron kinetic 
energy is not equally distributed in all directions arising from electrons in the well obtain extra kinetic 
energy in z-direction from the potential barriers. These are important findings which were overlooked in 
previous works. 1_ “ 


Second, the larger scattering angle 0 yields a larger value of For q> 1 , i.e., for the scattering 

processes in which electrons change their initial momentum direction, the value of l33 is almost independent 
of 0 and can be distinguished from various values of r by the second term of the right hand side of Eq. (20). 
The value of lQ 2D is independent of 0 and insensitive to the L z for q<i giving rise to a weak L, dependence 
or total electron scattering rates. It also appears that the value of I® 2 ® for q>l is rmre sensitive to L 
for q<1. z 


than 


Third, a dent occurs in each dotted curve which teco-.es smaller when r increases. The position of the 
dent moves towards the smaller q when r decreases. It should te noticed that the dent does not * appear in 
solid curves for Q2D structures. This is because the non-commutation relationship given by Eq. (1) employed 
for Q2D systems which gives rise to the momentum change in the z-direction does 
momentum change in the direction parallel to the quantum well interfaces. 


correlate to the 


The total normalized phonon emission rate (W k . k /* 0 ) in Q2D quantum wells versus nor-alized energy E k /huj 
Is displayed in Fig. 3 by two solid curves. The normalization factor «„ is given by Wo^rlo/h’q,*. For" GaAs 
W»-1.2k ps" 1 . For comparison, the calculated results for Q2D by Leburton" (dash curves) and for" 3D '(dotted 
curves) are also plotted in Fig. 3. The results of our approach clearly differ from the others 1 ** in two 


aspects. First, the Q2D phonon emission rates for various L z are 3 times smaller at the onset of phonon 
emission (E k /hw 0 -l) than previously obtained. 1 ' Second, the phonon emission rate is a slowly increasing 
function of Initial electron energy instead of being a decreasing function. The non-zero value of the phonon 
emission rate for E k -hu () stems from the nature of Q2D density of states. For thick quantum wells (L >5L 0 ) 


the phonon emission rate is about a factor of 2 smaller than the 3D rate for E k /hu 0 >2 resulting from the 
weakening of the electron-POP interaction in the z-direction for Q2D case and neglecting higher intra-subband 
and lnter-subband contributions. The Q2D emission rate will approach the 3D result when L z ->» if all the 
contributions mentioned above are included. The phonon emission rate is comparable to the 3D rate for very 
thin quantum wells arising from the competition between Q2D density of states and the weakening of electron- 
POP interaction strength in z-dlrectlon. 
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For completeness, it should be pointed out that the phonon absorption scattering rate can be treated in 
the same way as outlined in this paper. One should expect that the emitted phonon by Q2D electrons 
traveling in the direction perpendicular to the quantum well interface may be easily reabsorbed. This is 
because those phonons are not plane-wave-iike phor.on mode that is extended to the whole sample-. The 
phonons are localized inside the well the same as the localized electrons. 


In conclusion, the confinement effects on the phonon emission rate of the electron in Q2D structures 
through tne unscreened electron-POP interaction are shown by taxing into account the non-commutation 
relationship between the Hamiltonian of the system and the electron momentum operator in the z-direction. 
It is explicitly remonstrated how one can exactly treat the fuzzy momentum conservation in z-direction using 
neither "non-momentum-selection rule" nor "momentum conservation approximation" by considering the phonon 
wave vector q z as a r.v.. The phonon emission rate at the onset is found 3 times smaller than previously 
reported and increases with initial electron energy. 
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Fig. 2 The electron-PCP interaction strength for Q2D and 3D systems versus the parallel component of phonon 
wave vector q. The solid curves from 1 to » are for the 02D structures with well width 0.5 L 0 , L 0 , 3L 0 , and 
5l 0 accordingly. The main contribution of the integral I Q25 to the total phonon emission rate occurs at q-1, 
where lQ 2D presents a weak L z dependence. The dotted curves from 5 to 7 are for the 3D system with ratio 
r*10, 3, and 1.2, respectively. The scattering angle 8 for those three curves are set to be all zero. The 
dotted curve 8 is also for 3D system with r-1.2 and 8-ir/6. The larger scattering angle 8, the larger the 
value of I^ D in the range from 0.1 to 1 for q. 
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Optical Transitions and Enhanced Recombination Lifetimes in 
Quasi-zero Dimensional Electron System in CdS x Se,_ x 

Kai Shum, G. C. Tang, Mahesh R. Junnarkar, and R. R. Alfano 

Institute for Ultrafast Spectroscopy and Lasers 
Depts. of Physics and Electrical Engineering, City College of New York 
New York, New York 10031 


Abstract 

The recombination lifetimes for the radial and angular quantum number conserving 1S-*1S and 1P-»1P transi¬ 
tions from three-dimensionally confined electrons in CdS x Se,_ x were measured by time-resolved photolumines¬ 
cence. The assignment of the observed transitions is supported by calculations of eigen energy levels, 
squared matrix element ratio for those transitions, steady-state and picosecond spectroscopic studies. 

Introduction 

Carriers localized in semiconductor microstructures have attracted much attention both for fundamental 
reasons and for their potential applications. Semiconductor crystallites in a transparent insulating matrix 
can confine electrons in three dimensions because their diameters range from 30 X to 800 X.'“* The envelope 
wave functions in spherical coordinates [r, 0, $] and eigen energies of conduction electrons and valence 
holes localized in an infinite-spherical-well within the effective mass approximation are given by: 

|nLm > - C nL j L (xnL r > *Lm (6.*). 


and 


_ 2 Xn 2 L 

2m e ,h a 2 ’ 


( 2 ) 


respectively. The subscripts n,L, and m are effective radial, angular and magnetic quantum numbers, respec¬ 
tively. j and Y are the spherical Bessel and spherical Harmonic functions, respectively, a is the radius of 
crystallite and m e p is the effective mas3 of electron or isotropic hole mass. 5 The value of xnL *' or the 
lowest two states 'of either conduction electron or valence hole IS and IP are n and b.h9, respectively. S 
stands for L-0 and P for L*1. The allowed transitions which conserve angular and radial quantum numbers are 
1S-»1S, 1P-»1P, and higher transitions. The physical picture of the quasi-OD electron system described above 
has been experimentally verified by several groups. 1-5 Large and fast optical nonlinearities arising from 
photogenerated electron hole pairs have been also reported in quasi-OD electron systems."" 7 Pioneering sub¬ 
nanosecond time-resolved photoluminescence" and picosecond pump-probe experiments’ have been performed. It 
is important to study the dynamics of various excitations in quasi-OD electron systems in detail in order to 
understand the physical origins underlying these various novel optical properties. 


In this paper, we report on measurements of ultrafast recombination lifetimes for the radial and angular 
quantum number conserving 1 S-» 1 S and 1P-»1P transitions in quasi-OD electron systems in CdS x Se,_ x at . 3K 
using a streak camera detection system, and on the observation of optical transitions using photoluminescence 
spectroscopy between quantized levels (IS, IP) in the conduction and valence bands in quasi-zero dimensional 
electron system in CdS x Se,_ x . The time-resolved photoluminescence detected at various emission ener¬ 
gies allows us to unambiguously identify the 1S-1S and 1P-»1P transitions. The recombination lifetime 
of 1P-*1P transition was measured to be 3-5 times shorter than 1S-1S transition. The ultrafast decay of IP 
excitation may have practical importance for the construction of ultrafast reversible optical switches. The 
assignment of the observed S and P transitions were supported by calculations of eigen energy levels 
and squared matrix element ratios for these transitions, as well as from steady-state photolumlnes- 
cence spectra of the same samples used for the picosecond kinetic studies. 


Ex£e r i_m enta_l _M etho d 

The samples Investigated were four optical glass filters 2-6!, 2-59, 2-58 and 2-6h manufactured by Corn¬ 
ing years ago. The samples were mounted on a cold finger of a variable temperature helium optical cryostat 
(2K to 300K). The x value for each sample was accurately obtained from chemical analysis. 10 The picosecond 
experimental setup used in this research has been described in detail elsewhere. 11 A second harmonic (530 
nm) of a Nd-gla3s laser pulse of 8 ps duration was used to excite the samples on the front surface. The 
maximum optical energy Incident onto the front surface was about 90 pj. This corresponds to 10 1 * photons 
per single shot. The spot size was 8x10~* cm 2 . The exact photoexclted carrier density in quasi-OD electron 
system was difficult to estimate due to the complexity of the surface of the crystallites. Assuming all the 
photons were absorbed with 10pm, a value on the order of 10*" cm" 5 was estimated. The photoluminescence 
emitte : from the samples wa3 collected by a combination of lenses and imaged onto a 30-pm slit of a Hama¬ 
matsu streak camera system. Various narrow-band filters centered at different wavelengths were placed in 
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front of the streak camera to select different transition energies. Two color filters 3*67 were placed in 
front of the streak camera to block the light with wavelengths shorter than 5 *JOoJL The photolominescence 
profiles were corrected for the nonlinearity of streak rate and spectral response of the detector (S-20). 
For steady-state photoluminescence experiments, an Argon-ion laser, a double Spex spectrometer, 

S-20 photomultiplier and a lock-in amplifier were used. 


Picosecond Luminescence Studies 


The photoluminescence profiles in time domain were measured at different energies covering the comes 
ponding spectral range of each sample. The expression for number of photons detected at energy £ at time 
is given by: 

-t'saiE) _ 


N(E,t) = N„[e‘ 


.(E) 


], 


where E is the central energy of the i.arrow-band filter used; a d (E) and t r (E) are the decay and the rise 
times at energy E; and N 0 is the proportionality constant. The value of t r gives information about how fast 
electrons and holes scatter out from initially excited states to iess energetic states where they recombine 
while, gives a measure of the recombination lifetime of various excited states when recombination domi¬ 
nates the depopulation of photo-excited carriers. 



The time-resolved photoluminescence profiles obtained at various emission energies at H.2K for sample 2-58 
are shown in Fig. 1(a-e). The left curve in Fig. la shows the temporal profile of exciting laser pulse which 
reflects the time resolution of detection system (10 ps). The dotted curves in Fig. l are the fit to the 
data using Eq. (3) with a value of 12 ps for the risetime t r . The decay times of these time-resolved lumi¬ 
nescence profiles show only two distinct values. For emission energies ranging from 2.213 to 2.175 eV and 
from 2.1 to 1.967 eV the decay times are 29 and 100 ps, respectively. 


The decay time as a function of emitted photon energy E fo 
Fig. 2. For a comparison (Fig. 2e), the exciton lifetime 
CdS 0> 5 iSe 0 _», while the other two samples do not show tms f< 

reached In 2-61 and 2-59 samples). It is expected that t 


four quasi-OD electron systems are plotted in 
s exciton energy in quasi-ID system in bulk 
ure (the reason, the P levels have not been 
decay time should be nearly invariable over 


the energy range of emission if only the lowest confined state 13 for electrons and holes is occupied. The 
wide spectral range most likely reflects the crystallite size distribution and the fluctuations in the value 
of x from crystallite to crystallite. When the two lowest states 13 and IP are substantially occupied and 
the recombination lifetimes associated with these two states are considerably different, a steplike change of 
decay time will be observed in two distinct energy regions. The exact lifetime ratio in the two energy 
regions depends on the transition matrix elements. The above argument should in principle explain what we 
observed in Fig. 2. However, the energy dependence of idealized exciton lifetime in the bulk compound 
CdS„ SJ Se 0 ,, as displayed in Fig. 2e also exhibits a stepiike feature. The explanation given by Kash et al.* 
was "based" on the model suggested by Cohen and Sturge" where the exciton migration from a site with higher 
energy into another site with lower energy. This raises a question whether the exciton migration mechanism 
can apply to the present case. Since exciton migration must involve a transport over a relatively large dis¬ 
tance due to the nature of the exciton-phonon interaction," such migration may occur in relatively large 
crystallites. This apparently may explain the results observed in the sample 2-58 and 2-6H. In smaller 
crystallites the migration does not occur in correspondence with the samples 2—61. However, the exciton 
migration picture cannot explain the steplike energy dependence of lifetime as well as the ratio with the 
value of about 3.5 because the depth of a potential well 12 caused by compositional fluctuations are entirely 
random in nature. Furthermore, the exciton lifetime would be increased monotonically with a decrease of 
exciton energy as recently observed for CdTe-ZnTe interfaces. 1 ’ 


A two-state (IS, IP) model is introduced here which consistently interprets all the results shown in Fig. 
2 as well as other results observed by Cohen and Sturge." It is possible that the excitons are confined by 
the potential wells in the bulk alloy compound which have two lowest states designated by IS and IP. The 
exciton in higher-lying state IP can be scattered to tbe lower-lying state IS by emission of phonons instead 
of migrating from site to site. The two-state model can explain not only the results of Cohen and Sturge" 
but also our time-resolved dynamic results presented here as well as the data reported by Kash et al..* 


In order to further verify the two-states (IS, IP) model, we have calculated the eigen energy levels of 
the qua3l-0D electron systems and the recombination lifetime ratio for 1S-1S and 1P»1P transitions. The con¬ 
finement energy for the ground states is given by 


E >S ■ E ?S * El ?S 


(«) 


for the qua3l-0D electron systems, where Efg and E^g are the lowest confinement energies for electrons and 
holes, respectively. These energies were calculated from the bulk bandgaps'" at 300K and the measured peak 
energies of the fir3t derivative of room temperature reflectance (dR/di). 15 Using the electron effective 
mass and the isotropic hole effective mass given in Ref. R, the effective diameter 1 * of crystallite for each 
sample and eigen energies were calculated. These results are displayed in Table 1. The spin-orbit split 
states 1 " were neglected since they were not excited in our experiments. The emitted photon energies at 9.2K 
for 1S-MS and 1P«1P transitions were obtained by adding corresponding confinement energies to the values of 
bulk gap 17 at H.2K. The binding energies of excitons were not included. The calculated energy positions are 
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in good agreement with what we measured. These positions are located in Fig. 2 as arrows labeled by S and P 
for clarity. 

In Table 1, the values of lifetime for 1S-»1S transitions denoted by 1,3 were measured at the luminescence 
intensity peaks. A value of 85 ps for the 1S-»1S ground state recombination time for sample 2-59 is in exact 
agreement with the value obtained by Warnock and Awschalom* using pump-probe experiment. Since the energy 
per pulse in our experiment is about four orders of magnitude larger than they used, the Auger mechanism of 
carrier depopulation is then ruled out in these quasi-OD electron systems. Furthermore, we did not observe 
the change in lifetimes By varying the excitation intensity over a factor of 33. One concludes that the dep¬ 
opulation of photoexcited carriers must be dominated by radiative recombination path in the present study. 
The recombination lifetimes for 1P-* 1P transitions denoted by T,p are also entered in Table 1. Since the pre¬ 
sent experimental condition cannot detect the IP-OP transitions in samples 2-61 and 2-59 the numbers inside 
the circles are the anticipated lifetimes of 1P-»1P transitions assuming the same ratio of T,g/fip for the 
samples 2-58 and 2-69. Note that a value of 18 ps has been also reported by Warnock and Qwschalom’ for 
2-61 although they did not identify this to be the lifetime of 1P-*1P transition. It should be emphasised 
that strongly enhanced transition probability which leads to shortened lifetimes of various excitations in 
the quasi-zero dimensional electron system is r.ot a surprising result. The recombination enhancement due to 
carrier confinement in quantum well structure (quasi-two dimensional electron system) was clearly demon¬ 
strated by Gobel, Jung, Kuhl, and Ploog.'* These authors'* attributed their measured recombination lifetime 
to be the spontaneous carrier lifetime in quasi-two dimensional the electron system. The measured diameter 
dependence of recombination lifetime supports their arguments.'* 

A further test to support the two-states picture is whether a reasonable prediction can be made to locate 
the energy positions of exciton localized in IS and IP 3 tates for the quasi-OD electron system in the bulk 
compound CdS 0t s jSe 0 . * 7 studied by Kash et al..* When electrons are confined in a one-dimensional semiconduc¬ 
tor quantum well with thickness L z , the electron-hole recombination time for ground states may be estimated 
to be proportional to the L z since the wave function is "squeezed" by the boundaries' The evidence for 
this has already been provided by the measurements of Gobel et al..‘* Extending the above concept, the 
effective diameter (d) for the spherical potential well may be estimated using a relation T,g * c d*. A 
value of about 250 A for the diameter of the potential well in the bulk was estimated using the lifetime for 
energies less than 2.09 eV assigned by us to be the recombination lifetime for 1S-1S transition. Using this 
value the confinement eigen energies for exciton (electron and hole Dound By their mutual Coulomb potential) 
were calculated. These results are also entered in TaDle 1 under Duik and displayed in Fig. 2e. The experi¬ 
mental results of Kash et al.* are in excellent agreement with our two-states model. Furthermore, when the 
temperature of their* sample was raised above 3oK (3 new) the steplixe feature of lifetime in two distinct 
energy regions disappeared. Since exciton binding energy in this material is about 20 meV, the exciton disso¬ 
ciation is not responsible for this disappearance. However, the energy difference between ground state IS 
and first excited 3 tate IP for hole is just 3 meV (Table '); therefore, the strong hole scattering between IS 
and IP 3 tates assisted by thermal energy is responsible for the absence of the two distinct regions above 
36 k. This temperature dependence of exciton lifetime further supports our model. 

The lifetime ratio for 1S-• 1S and IP-lP can be predicted following the approach given by Casey and Panish” 
and Efros and Efros. 20 The matrix element for 3 transition from localized state in conduction band to the 
localized state in valence can be expressed as: 

■“* * -b H env> ( 5 ) 

where is the average matrix element for the 31och states for bands in absence of new eigen states due to 
the confinement and M env is the envelope part of the matrix element which is given by: 

^env *■ <1L | 1L>. (6) 

Since IP state is a threefold degeneracy state with magnetic quantum number -1, 0, and 1 , a value of 5 for 
the ratio of |M env |i^ ( p and |M env li’s-^s can be readily calculated. This value should be compared with 
experimentally measured value of 3.5 for t, 3 / 3 , 0 . The agreement is reasonable in view of assumptions made 
in the theoretical approach and the uncertainty in data fitting. 

St e a dy_j s t a te_ _P h 0 to 1 urn i_n_e_s_cence_ _St u di_e s 

Following work of Warnock and Awschalom’ steady-state photoluminescence studies was used to confirm our 
findings. In Fig. 3a, photoluminescence spectra of sample 2-6] at various excitation levels taken at room 
temperature is plotted. The excitation wavelength was the 988 nm line of the Argon-ion laser. The excita¬ 
tion levels were 2.5 and 10 times higher for the broken curve and dot-dashed curve than the solid curve 
respectively. The scale of luminescence intensity was 2 times and 10 times larger for the broken curve and 
dot-dashed curve than the solid curve, respectively. The broken vertical line indicates the peak position of 
dR/dA. The most salient feature we would like to stress is the appearance of the two-peak structure at the 
high energy side of the peak position of the first derivative of room temperature reflectance spectrum. 
Based on results of picosecond studies listed in Table 1 we attribute these two peaks to the 1S-1S and 1P-.1P 
optical transitions* as indicated by S and P in the figure. The energy separation between P and S is about 
220 meV. This rules out the possibility of transition from IS electron state to IS spin-split-hole 3 tate, 
since the energy separation of spin-orbit splitting is about 350 meV for Se mole composition of 0.73. Two 
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other features should also be mentioned. The excitation dependence of luminescence intensity did not show 
the saturation of available electronic states up to the maximum excitation level lCGW/'cm 2 : *nd the lumines¬ 
cence peak shifts to the lower energy when the excitation level increases reflecting the renormalization 
process in quasi-zero dimensional electron system. 

The photoluminescence spectrum of 2-61 measured at ^. 3K using the “57.9 r.m line of the Argon-ion laser as 
excitation wavelength is shown in Fig. lb. As can be seen from the spectrum the position of 1S-1S peak 
snifts to the higher energy and its FWHM (58 meV) decreases by a factor of 2 in comparison with room tempet— 
ature data. The latter indicates the broadening of 1S-1S transition at room temperature is not entirely 
dominated by size fluctuation of the microcrystallites under the photoexcitation. The line broadening 
mechanisms may apply for the quasi-zero electron system. Vie assume the broadening of the IS peak is domi¬ 
nated by the fluctuation of the diameter of the microcrystallites at 9.3 K. The inset of Fig. lb shows a 
spectrum fit by assuming the Gaussian shape 21 probability density function for the random variable d. 21 The 
fit gives a mean value of 79 X for the effective diameter and a variance a 2 - 25 « which characterizes the 
size fluctuation of microcrystallites. The structure beside the main IS peak is from the 1P-1P transition 
since its energy separation with the IP peak is about 210 meV. The broad peak centered at 507 n.m may arise 
from the transition from the split orbit band. 

In Fig. 9, the photoluminescence soectra of sample 2-59 measured both at room temperature (a) 3nd 9.JX (b) 
are displayed. The most important feature is the appearance of both 1S-* 1S and lP-*ip transitions. The 
energy separation between IS and IP peaks is 219 meV. The calibrated value for the confinement energy is 
213 meV. The corresponding peak positions are shifted to lower energy in comparison with 2-61 due to the 
larger diameter of crystallites and smaller mole percent of the composition of sulfur. 

The broad peak centered at 500 nm in Fig. 9b at the 9.3K spectrum appears in the 9.3K spectrum of the 
2-61 sample. The origin of this emission peak at 500 nm is unlikely from the ID-ID transition. This emis¬ 
sion may arise from the optical transition from the IS electron state to the IS spin-orbit split hole state. 
The energy separations between the weak peak and the main IS peak are about 370 meV and -50 neV for 2-61 
and 2-59, respectively. These energies are consistent with the measured selenium-composition-percentage dep¬ 
endence of spin-orbit splitting. 1 * 

We have also measured the room temperature spectra of 2-58 and 2-69 samples. The peaws for 2-58 and 
2-69 are at 1.9537 eV and 1.9119 eV, respectively. These peaks were broad so that the IP structures were 
not clearly resolved. 

Sumrnar y 

In summary, we have reported on the measurements of recombination lifetimes of 15-IS and 1P-JP transi¬ 
tions in quasi-OD electron systems by time resolved photoluminescence measurements. The assignment of these 
observed transitions are supported by calculations of eigen energy levels and the corresponding matrix ele¬ 
ment ratio as well as the direct observation of optical transitions between quantized energy levels (IS, IP) 
In the conduction and valence bands by steady-state photoluminescence measurements. 
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Table 1: The measured energy peak of first derivative of reflection at 300K, bulk energy gaps at k.2K and 
JOOK, the calculated confinement energies, emission energies for IS -► IS and IP -► IP transitions and effec¬ 
tive diameters (d), the measured lifetimes for IS -» IS and IP ■* IP transitions as well as their ratio are 
lisplayed. The t are anticipated values assuming a same ratio T, 3 /T,p for the samples 2-58 and 2 - 6 ^ . 
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Fig. 1 Time-reaolved photolumlnescence profiles detected at different energies (a to e) for the sample 2-58 
at n.2K. The dotted curves are the theoretical fit with 12 ps for rlsetlme and decay times are indicated on 
the corresponding profiles. The left curve of (a) is the temporal response profile to the laser pulse used 
to excite the samples. 
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Fig. 3(a) Photolumlnescence spectra of sample 2-61 at various excitation levels taken at room temperature 
using «88 nm line of the Argon-Ion laser. The excitation levels (Intensity scale) were 2.5 (2) and 10 (10) 
times higher for the broken curve and dot-dashed curve than the solid curve, respectively. The broken verti¬ 
cal line Indicates the peak position of the first derivative of reflectance at room temperature. S and P 
stand for the 1S-*1S and 1P-»1P transitions, respectively. (b) Photolumlnescence spectra of 2-61 at A.3 k. 
The Inset shows that the IS peak broadening Is dominated by the diameter fluctuation of microcrystallites. 
The solid dots are calculated by assuming a Gaussian shape of the probability density function for the random 
variable d with a variance of 25 A. 
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Fig. 4(a) Photoluminescence spectra o' 2-59 at room temperature excited by 488 nm line of the Argon-ion 
laser. The excitation level for the broken curve is half of the solid curve while the intensity scale is the 
same for the two curves. S and P stand for 1SS and 1P»1P transition, respectively. The broken solid line 
indicates the peak position of the spectrum dR/d (at room temperature). (b) Photoluminescence spectrum of 
2-59 at 4.3 K. 
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The salient points of our previous paper [Phys. Rev. B 33, 7259 (1986)] are reinforced in this 
Reply to the preceding Comment by Miller. 


The essential points of our work 1 were to demonstrate 
what experimental data from optical transitions should 
be used and how sensitive they are to the Q value of band 
offset. Contrary to Miller’s comment. 2 Dingle et al. 1 did 
not emphasize how sensitive the Q value depended upon 
the value of the energy separation, A E(L : ), between the 
n =1 heavy- and light-hole subbands. We found 1 that 
the values of A E are most sensitive to the band offset for 
L, ranging from 15 to 80 A. Based on this key finding a 
currently fashionable choice of £),, =0.40 proposed by 
Miller et al.* was ruled out using Dingle’s connection 
rule and the then available experimental data in this 
range of L.. Another key point made in our paper 1 was 
that systematic measurements should be performed in 
the sensitive zone of well width (15 A to 80 A) to pre¬ 
cisely determine the Q value. These essential points are 
still overlooked in depth by Miller's comment 2 to our pa¬ 
per. 1 

It should be pointed out here that the extent of data 
used in our paper 1 was not the main issue of the paper 
and more data could not alter the essential spirit of the 
paper. Most recently, Miller et al. have performed new 
measurements providing more data in the sensitive zone 
which will be discussed below in this Comment. 

We have recognized that there has been a large body 
of data pointing to a larger Q„~0.40, which was the 
mainstream of thought in this field. However, most of 
the experimental probes carry their own source of uncer¬ 
tainties. 6 The electrical measurements are often plagued 
by residual doping whereas the intersubband optical 
transitions may prove to be too weakly dependent upon 
the band offset. The optical method described in our pa¬ 
per 1 provides a sensitive optical test of band offsets in 
GaAs/Al,GaAS|_, structures. It should be mentioned 
that there are some most recent determinations of Q 
value for GaAs/AI^Ga^^As which do not yield 
g„=0.40 but point to values of 0.31 from Raman 
scattering data 7 and 0.23 ±0.07 from electrolyte 
electroreflectance study.* 

As noticed by Miller, 2 one of the data in the original 
paper 1 to support lower Q v was obtained in a circular 
fashion. However, the other data quoted by us do not 
suffer from this problem. 

It is well known that the mismatch of effective mass in 
well and barrier is responsible for the appearance of 


various connection rules'' 9 ' 10 for the envelope wave func¬ 
tion and its derivative. In general, using a different con¬ 
nection rule will result in different eigen energy levels 
and hence differences in the spacing between the levels. 
Therefore, the Q value determined by A E versus L. de¬ 
pends much on what connection rule is employed. This 
is where some of the problems come from. 

In the spirit of our past paper, the light- and heavy- 
hole energy splitting is plotted in Fig. 1 as function of 
well thickness L z using various connection rules. 3,911 
The Q value is adjusted in such a wav the A£s for 
L z > 40 A calculated 12 using different connection rules 
coincide with each other within —2 meV. D, B, and ,Y 
in the Fig. 1 refer to the Dingle’s connection rule, 
Bastard's connection rule, and new connection rule. 11 
The material parameters used for the calculation of the 



L z (A) 


FIG. 1. Calculated light- and heavy-hole energy splitting as 
function of L,. D. B, and .V refer to Dingle’s, Bastard's, and 
new connection rules. The values of Q„ are indicated on the 
corresponding curves. The masses used for the dotted curve is 
Miller’s masses. The masses used for the other curves are the 
conventional masses given in the text. The plus and triangle 
data are from Refs. 5 and 15, respectively. 
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curves are jc= 0.37, A£T,. = 1247jc0 (r (meV) (the value of 
Q,, is indicated next to the corresponding variable in the 
Fig. 1), masses' 1,14 of heavy hole (light hole) for GaAs 
and AlAs are 0.403m,, (0.087m (i ) and 0.487w fi 

(0.208m o ), respectively; except for the curve where the 
heavy-hole mass (0.34m 0 ) and light-hole mass (0.094m„) 
proposed by Miller et al.* were used. The plus and tri¬ 
angle data points in Fig. 1 are from Refs. 5 and 15, re¬ 
spectively. 

Figure 1 shows two important features. First, as ex¬ 
pected the energy splitting A£ is sensitive to the connec¬ 
tion rule in the sensitive zone of L z , especially for 
L z <40 A. This features makes it possible to experimen¬ 
tally demonstrate which connection rule is appropriate 
by systematically measuring A£ in the sensitive zone. 
Based on the data given by Miller et al., s,l> Dingle’s 
connection rule which was used by us to show how to 
determine Q value in the sensitive zone of L t seems to be 
not appropriate. However, the method itself described 
in our paper 1 is still very useful. Second, for the given 
connection rule the splitting A£ is very sensitive to the 
effective masses. In order to fit simultaneously the data 
given by Miller et al. 5,15 it seems to be necessary to use 
the Bastard’s connection rule and the mass parameter set 
< m HM =0.34/n 0 , m LH = 0.094m 0 ) and Q ,. = 0.40 pro¬ 
posed by Miller et al.* The use of new connection rule 11 
with £>„=0.23 seems to fit the data for L z > 30 A. It 
should be pointed out here that in our calculation the 
binding-energy difference between the light hole and the 
heavy hole is taken to be a constant ( —0.5 meV). This 
may be a good approximation for L z > 40 A. In order to 
comparethe measured A £ with the calculated A £ for 
£,<40 A an exact analysis taking the proper exciton 
binding energy into account must be performed. 16 

In Fig. 2, the energy separation A£ 12 between first 
and second conduction subbands is plotted as function of 
L t using three different connection rules by including 17 
the energy-dependent mass m„.(£„) (Ref. 18) as well as 
m fc (£„ ), n = 1,2, for well material GaAs and barrier ma¬ 
terial Al 0 jGa 0 7 As, respectively, and compared A£, : 
with the published experimental data by West and 
Eglash. l<> This work is not affected by the valence-band 
complexity, the exciton binding energy, and can provide 
an independent test of Q value. 

The m b (E„ ) is obtained by 


a 2 

F 2 

ZL + 

fi 2 k 1 E t 

1/2 £ 

a 2 

* 2 k 2 

a* 2 

4 + 

2m CD 

2 

dk 2 

2tn b (E „) 


( 1 ) 

The expression inside of the large brackets on the left- 
hand side of Eq. (1) is the energy-momentum relation¬ 
ship within semiconductor band gap. 20 After 
differentiation the m b (E„) reduced to 

m b {E„ )=m CB ( 1 +4V, ) V2 , n = l,2, ... , (2) 

where m C8 is electron effective mass at conduction band 
edge of AI 0 J Ga 07 As, X„ = I ff 2 fc„ 2 „ /2m c „£ ( |, and 
n = l,2 for first and second subbands. Note that k 2 is 
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FIG. 2. Calculated energy separation A£ i; between first and 
second conduction subbands as a function of well width L. us¬ 
ing Dingle's ID), Bastard's IB), and new (,V) connection rules. 
The curves are calculated using A£ r =(? f A£, (A£, = 1247* 
meV, jc = 0.30). The values of Q c arc indicated on the corre¬ 
sponding curves. The data of West and Eglash I Ref. |9> are in¬ 
dicated as two plus signs. The inset shows schematically the 
well and the energy-momentum dispersion curves. The elec¬ 
tron effective masses at band edges for GaAs and Al„ ,Ga 0 ,As 
are also indicated in the inset. 

negative and k b is imaginary reflecting the electron wave 
function within the gap of Al 03 Ga 07 As is exponentially 
damped. The dispersion relations dependent upon the 
real and imaginary wave vectors in the well and the bar¬ 
rier are schematically shown in the inset of Fig. 2. The 
connection rule and the Q value used are indicated on 
the corresponding curves. The value of .t is equal to 
0.30 in order to make a comparison with the data A£, 2 
of West and Eglash. 19 

Several features are very apparent in Fig 2. First, for 
the given connection rule A£ l2 is very sensitive to the Q 
value in the sensitive zone of L. giving rise to a possibili¬ 
ty of extracting (he value of A£ c and Q value. Second, 
the value of A£ l2 calculated using the new connection 
rule approaches the values of A£ ]2 at —45 A and at 
L z > 100 A calculated using Bastard's and Dingle's con¬ 
nection rules, respectively. Third, the two experimental 
points support neither Bastard's nor Dingle’s connection 
rules, but seems to be more favorable to the new connec¬ 
tion rule with Q c =0.77. 

In conclusion, we have shown different connection 
rules and hole masses yield different Q values. Further¬ 
more, the appropriate connection rule can be discrim¬ 
inated by systematically measuring A £ and A£, 2 in the 
sensitive zone of L z . Dingle’s 85-15 rule is based on his 
connection rule while Miller’s b0-40 rule is based on 
Miller's effective masses and Bastard's connection rule. 
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The essential.information described in our paper 1 is not 
compromised by the comments raised by Miller but rein¬ 
forced here. We stress, that more systematic photo- 
luminescence r PLE, and inter-conduction-subband tran¬ 
sition data must be taken in the sensitive zone (15 to 80 


A) to determine the proper connection rule and an accu¬ 
rate Q value. 

This work was supported by Air Force of Scientific 
Research under Grant No. AFOSR-86-0031. 
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Gallium arsenide photoluminescence under picosecond-laser-driven shock 
compression 

X. Z. Lu. R. Garuthara. S. Lee. andR. R. Alfano 

Institute for Ultrafast Spectroscopy and Lasers. The City College of Mew York, Mew York. Mew I ork 10031 
(Received 8 September 1987; accepted for publication 10 November 1987) 

A pump-and-probe technique was used to investigate shock effects on the photoluminescence 
spectra ( —833 nm) at T = 80 K due to the direct transition £ 0 from the T 6 conduction band 
to the r„ fourfold degenerate top valence band in GaAs. Under the shock loading condition, 
the photoluminescence peak was observed to blue shift and split into two components, 
corresponding to the transitions from the r„ conduction band to the valence heavy- and light- 
hole subbands, because of symmetry breaking by the uniaxial shock compression along the 
[001 ] direction. From the blue shift of the photoluminescence peaks, we deduced our 
picosecond-laser-driven shock pressure of ~ 10 kbar. 


Significant changes in the optical and transport proper¬ 
ties of semiconductors occur with the application of a uniax¬ 
ial stress due to changes in symmetry and lattice param¬ 
eters. 1-1 Knowledge of stress effects on band structure is 
important for operation of GaAs in high-speed ultrasmall 
devices and switches. In this letter we report on a new obser¬ 
vation of the photoluminescence spectra from gallium arsen¬ 
ide (GaAs) under picosecond-laser-driven shock wave load¬ 
ing. The shock waves were generated by focusing intense 
picosecond laser pulses onto an aluminum (Al) foil 5 '" at¬ 
tached to the sample. 

GaAs has a zinc-blende structure 1 ' of the space group 


Tl. The lowest conduction band is s-like and twofold degen¬ 
erate, having a T 6 symmetry at the Brillouin zone center. 
The top valence band is a sixfold degenerate/>-like band that 
splits into a fourfold P }/2 multiplet ( J = 3/2, Mj — + 3/2, 
± 1/2) with a r„ symmetry and P in doublet (J = 1/2, 
Mj = +1/2) with a T 7 symmetry due to the spin-orbit in¬ 
teraction. The band structure is shown in Fig. 1 (the left- 
hand side). The lowest energy band-to-band direct transi¬ 
tion (labeled £ 0 ) occurs at the center of the Brillouin zone. 
The total Hamiltonian matrix for the valence band in the 
presence of a uniaxial compression P along the [001 ] direc¬ 
tion can be written’' 


I 


H 



0 

- SE h + J SE s 
2~' n 6E s . 



( 1 ) 





1 

% 


£ 

i 


i 


where SE„ = — a (5,, + 2S U )P, SE H = — a'(5 n 
+ 2S lz )P, SE s = -2b(S u ~S i2 )P, SE s 
— — 2A'(5 n —5,2)/*. and A 0 is the spin-orbit splitting, 5, i 
and S, 2 are components of the elastic compliance constants, 
a and a' are the deformation potentials of the T g and T 7 
valence bands for hydrostatic effects, and b and b' are the 
deformation potentials of the T 8 and T, bands for shear ef¬ 
fects. If the pressure-induced shift and splitting are much 
smaller than A„, the eigenvalues of Eq. (1) can be written 


<S£„, = -<5£„ - (<5£ s , 

(2a) 

SE„ = SE H + 1<5£ S + (<5£v> , 

(2b) 

2A 0 


*£-. " Ao 6E h {SEs ) \ 

(2c) 

2A 0 



Including the effect of the conduction-band shift, the pres¬ 
sure dependence of the energy gap changes is given by 

6(E, - £„ ) = A£„ + J 6E S , (3«) 


(SE k ) 2 

S(E C - £„ ) = A£„ — J SE s -—, (3b) 

/cr \2 

S(E C — £„.) = A 0 + A£ w + ■ _ s , (3c) 

2A 0 

where A£ H = — (c, + j)(5,, + 25 u )£. A£„• 

= — (c, +a') (5,i + 2S U )P, andc, is the deformation po¬ 
tential of the conduction band. For light polarized perpen¬ 
dicularly to the stress axis c, the transitions between the r- 
like T 6 conduction band and all three valence bands, v,,v 2 , 
and v y , are allowed. 9 From Eqs. (3a) and (3b), we deduce 
that under shock loading condition, the photoluminescence 
peak corresponding to the E 0 transition will be subject to a 
blue shift because of the hydrostatic component of the pres¬ 
sure [represented by the term A E H in Eqs. (3a) and (3b) ]. 
It will also split into two components corresponding to the 
transitions £„( l) and£ 0 (2 ; (see Fig. 1, the right-hand side) 
because of the reduced symmetry by the uniaxial nature of 
the shock pressure [ represented by the term 8E S in Eqs. (3a) 
and (3b) J. Since the heavy-hole band in has a lower energy 
under uniaxial compression, 10 the corresponding transition 
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UNSHOCKED SHOCKED 

FIG. 1. Band structure of GaAs at T = 80 K for unshocked (left) and 
shocked (right) crystals. 


£ 0 (2) will have a larger blue shift. In addition, the blue shift 
of the emission peak of the electron-light-hole recombina¬ 
tion [the transition £„( 1) 1 will increase sublinearly with the 
increasing pressure because of the pressure-induced cou¬ 
pling between the sets of Mj ~ + 1 bands (v, and v,). 

A schematic diagram of the experimental geometry is 
shown in Fig. 2. The experimental setup was described in 
Ref. 11. A Quantel NdrYAG laser of 30 ps pulse width and 
25 mJ pulse energy at l .064 gm was utilized as the pump 
beam to generate shock waves in the aluminum foil. The 
shock waves were launched into the GaAs sample located 
behind the aluminum foil. The probe beam of 27 ps at 532 
nm was delayed through a white cell approximately by 46 ns 
relative to the arrival of the pump beam at the aluminum foil 
to ensure that the probe pulses reach the sample surface area 
right after the arrival of the shock front. This was necessary 
to observe the maximum effect of shock pressure. 

The GaAs sample was n type with a carrier concentra¬ 
tion of ~ 10 ,K cm~\ The thickness of the sample was about 
50/tm with the cubic [001J axis perpendicular to the sample 
surface. The aluminum foil had a thickness of 20 gm. The 
sample and aluminum foil were loaded between two quartz 
plates and placed in the center of a liquid-nitrogen Dewar. 



FIG- 2. Schematic diagram of the experimental geometry 


APIEZON N-grease was added to the interface of the sam¬ 
ple and aluminum foil to ensure efficient mechanical cou¬ 
pling at T — 80 K. The pump beam was focused by a 15-cm 
focal length lens to a 500-/im-diam spot on the aluminum 
foil, which ensured planarity of the shock wave and a peak 
power density of — 4.2X 10" W/cm : . The probe pulse was 
focused to a 400-/im-diam spot by a 15-cm focal length lens 
on the sample, giving a peak power density of — 1.6 X 10 7 W/ 
cm : . 

The emitted light was dispersed by a 1/4 m spectrom¬ 
eter, delected by a silicon-intensified target coupled to an 
optical multichannel analyzer and stored in a PDP 11/23 + 
host computer. The resolution of the detection system was 
about 2 nm. All spectra were recorded at T — 80 K by single 
shots of the laser. 

The observed photoluminescence spectrum of GaAs 
without shock loading is dominated by a single peak due to 
the band-to-band recombination at the lowest energy gap £,„ 
as shown in Fig. 3 (left). No emission lines due to impurity 
transitions were observed. The effect of shock waves was 
most pronounced for the delay time between the Dump and 
probe beam, r = 40 —60 ns. At r< 35 ns, no shock effect was 
observed. A typical spectrum showing the shock effects on 
the photoluminescence from GaAs is displayed in Fig. 3 for 
r = 46 ns and probe intensity 1.6 X 10 W/cm\ Under shock 
loading, the emission peak was split into two, with blue shifts 
of 20 and 85 meV. respectively. 

The lowest energy band gap in GaAs in the absence of 
the shock pressure can be written 1 -' 

E,lT) = 1.530eV - (50* 10 * eV K ')T. (4) 

The peak energy of the emission due to the band-to-band 
transition £„ at T — 80 K is E p — E„ (80 K) = 1.490 eV, 
i.e., A. r = 832 nm, which agrees with our observed un¬ 
shocked emission line at 833 nm ( Fig. 3) within our experi¬ 
mental accuracy. Since the energy gap E, — £, expands 
linearly with the increasing pressure P, the blue shift of the 
transition energy £„(2) was chosen to calibrate the shock 



FIG .V Unshocked (left) and shocked (right) photoluminescence spectra 
at a delav time r * 46 ns and probe intensity I h « 10' W/cm : The shocked 
spectrum consists of t*o peaks, corresponding to the transitions £,,( I) and 
£„( 2 ). 
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pressure. From Eq. (3a) the pressure coefficient of this tran¬ 
sition energy is given by 

dE„( 2) 

- dp = - (c, +fl)(S M + 25, 2 ) -b(S u -S l2 ). 

(5) 

With the reported values, 13 c,+a= —13.0 eV, 
6= -1.66 eV, and 14 S,, = 1.16X 10“ 3 kbar'\ 5,, 
= — 3.67x 10 -4 kbar - ', we obtain dE t) (2)/dP = 8.5 
meV/kbar. Therefore, the experimentally measured blue 
shift of 85 meV (Fig. 3) for the electron-heavy-hole recom¬ 
bination emission line gives a magnitude of the laser-driven 
shock pressure in GaAs of ~ 10 kbar. 

From Eq. (3b) and the reported values, 4 b' = — 2.47 
eV and A () = 0.34 eV, we obtain S(E C - E , : ) = 21.5 meV 
for a pressure of 10 kbar, which consistently agrees with our 
observed blue shift of ~ 20 meV under the same shock com¬ 
pression condition. Therefore, the shocked system is in qua¬ 
siequilibrium state. 

In conclusion, the picosecond-laser-drive shock waves 
offer a new technique to investigate properties of materials at 
very high pressures with a table-top laser facility. Using this 
technique, we have observed blue shift and line splitting in 
the photoluminescence spectrum of GaAs due to band-gap 
expansion and symmetry breaking by the uniaxial shock 
compression along the [001 ] direction. The pressure range 
of the shock waves generated in our experiment by a 25-mJ 
picosecond laser pulse is — 10 kbar. 


This work is supported by Air Force Office of Scientific 
Research 86-0031 and Office of Naval Research. 
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Shock-wave-induced collision broadening of the photoluminescence spectra 
in GaSe 
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New York 10031 

(Received 13 August 1987; accepted for publication 28 September 1987) 

Significant spectral broadening of the photoluminescence in GaSe under the picosecond-laser- 
driven shock pressure has been observed for the first time. The broadening of the spontaneous 
emission was found to be proportional to the shock pressure and attributed to a shock-wave- 
induced exciton collision mechanism due to the directional motion of particles in the shocked 
region. 


There are significant differences' - * between the ob¬ 
served physical and chemical properties of condensed mat¬ 
ter which is shocked and hydrostatically compressed at the 
same pressure and temperature. Production of the picose- 
cond-laser-driven shock waves offers a new powerful tech¬ 
nique to study the properties of semiconductor circuits and 
devices at very high pressures with a table-top laser facility. 
Recently, we reported 5 on the spontaneous and stimulated 
emission from GaSe at different excitation intensities and 
attributed both the spontaneous and stimulated emission to 
the same origin—exciton-exciton scattering process. In this 
letter, photoluminescence spectral line broadening observed 
under picosecond-laser-driven shock loading in GaSe is re¬ 
ported for the first time which is attributed to enhanced exci- 
tonic collision processes associated with the shock wave. 

A pump-and-probe technique was used to observe the 
shock effect on the photoluminescence emission in GaSe. 
The GaSe sample of 50/rm thickness with the c axis perpen¬ 
dicular to the layers was attached to an aluminum foil of 
thickness 20 ptm. A Quantel Nd:YAG laser of 30 ps pulse 
width anH 35 mj pulse energy at 1.064/im was focused to a 
450-fti.i-diam spot on the A1 foil surface to generate shock 
waves which propagated through the A1 foil into the GaSe 
sample. A 532-nm probe beam of 27 ps pulse width was fo¬ 
cused to a 350-/im-diam spot onto the GaSe surface. The 
photoluminescence from GaSe was dispersed by a 1/4-m 
spectrometer and detected by a silicon-intensified target 
coupled to an optical multichannel analyzer OMA III, and 
stored in a PDP1 1/23 4- computer. 

Typical spontaneous emission spectra from GaSe with 
and without shock loading are shown in Fig. 1. The salient 
features displayed in Fig. 1 are red shift of 28 nm and broad¬ 
ening of the shocked spectrum relative to the unshocked one. 

The emission due to exciton-exciton scattering can be 
expressed as (£*,£* ) — (hv,e — h). The spontaneous 
emission spectrum due to the exciton-exciton scattering pro¬ 
cess is given by* 



(£ — E x ) 2 + (na/e)E\ kT, 



K 

(i + (£/£,)]* 



(l) 


where p(E)dE is the number of photon modes in the crystal 
between energies E and £ + dE, wherep(£) cc £ 2 , £, is the 
energy of the Is exciton at K = 0, E b is the binding energy of 
excitons, E x = E t — E b ,T x is exciton temperature, a and e 
are the polarizability of excitons and dielectric constant of 
the material, respectively. The peak energy is located at 
E p = E g — 2 E b . Figure 2 shows a comparison between the 
observed spontaneous emission spectrum (solid dots) and 
the calculated one (solid line) from Eq. (1) with E x = 2.0 
eV, E b = 0.020 eV, T x = 500 K, and the coupling coeffi¬ 
cient na/e = 3.3 X 10~ 3 . 

It is a well known fact that the direct energy band gap of 
GaSe shrinks under pressure causing the red shift of the pho¬ 
toconductivity, absorption, and luminescence spectra. 7 '’ 
With the reported static values of 10 3E t /3P — —6.2 

meV(kbar) ~ 1 and 3E b /3P = - 0.6 meV(kbar) ~ 1 we ob¬ 
tain the (hydrostatic) pressure coefficient of the peak energy 
3E p /3P = — 5.0 meV(kbar)"'. There is evidence 11 show¬ 
ing that the stress anisotropy of the shock compression could 
be neglected in many cases. The red shift of 28 nm (84 meV) 
displayed in Fig. 1 deduces a shock pressure of ~ 17 kbar. To 
see if this approximation is reasonable, we also performed an 
independent measurement of the shock pressure in the same 
experimental setup by an jc-cut quartz transducer tech¬ 
nique. 12 The pressure of — 15 kbar was obtained which con¬ 
firmed our approximation. 

Significant line broadening of the spontaneous emission 
spectra under the shock loading was clearly and repeatedly 
observed. The magnitude of broadening was about 24 nm for 
a shock pressure of 17 kbar, as displayed in Fig. 1. It was 
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FIG. I. Observed unshocked and shocked spontaneous emission sp e c t r a 
from GaSe at room temperature and excitation intensity 17 MW/cm 2 . 
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FIG. 2. Calculated spontaneous emission spectrum from the exciton-exci- 
ton scattering model of Eq. (1) with E, = 2.0 eV, E„ = 0 020 eV, na/ 
f = 3.3 X 10 and T, = 500 K. The solid dots represent the experimental 
unshocked spontaneous spectrum. 

found that the peak amplitude decreased under shock load¬ 
ing such that the integrated spontaneous emission intensity 
(area) approximately remains the same as that of the un¬ 
shocked spectrum. 

It was reported 13 that although the electronic (or exci- 
tonic) temperature in the probed region (~ 434 K) was con¬ 
siderably higher than the room temperature due to heating 
effect by the probe beam, no significant temperature rise due 
to the shock wave loading was detected. The Doppler 
linewidth at temperature T x due to the thermal motion of 
the emitting excitons is given by 14 



where A p is the peak wavelength of the emission and M x is 
the exciton mass. If this mechanism was responsible for the 
observed line broadening (of a factor of 2, approximately), 
the temperature would increase by a factor of 4 due to the 
shock loading, i.e., the exitonic temperature would exceed 
1700 K in the shocked region which is definitely unrealistic. 
Also, a Doppler broadened line should have a Gaussian pro¬ 
file 14 which disagrees with our observed Lorentzian line 
shape. Therefore, the possibility of temperature effect caus¬ 
ing the linewidth increase due to the shock loading can be 
safely excluded. 

The probed region in our experiment formed a disk of 
350/im in diameter and 5/zm in thickness (the penetration 
depth of the 532-nm probe beam in GaSe). Since the shock 
pressure is proportional to the pump beam fluence, 12 trans¬ 
verse inhomogeneity of the shock pressure over the disk area 
(350/ 2 m in diameter) originated from the Gaussian profile 
of the pump beam (450/im in diameter) intensity should be 
less than 20%. There was also a longitudinal inhomogeneity 
of the shock pressure in the direction of shock wave propaga¬ 
tion due to pressure decay in space. In our picosecond laser- 
driven shock wave case, the shock decay time was mea¬ 
sured 13 to be —70 ns, corresponding to a decay length of 
— 140/im. Thus, the longitudinal variation of pressure in the 
shocked region was less than 0.5 kbar. However, the ob¬ 
served line broadening corresponded to a pressure variation 
from zero to twice the mean value of the shock pressure as 
deduced from the observed red shift. Therefore, neither 
transverse nor longitudinal inhomogeneity of the shock 
pressure can explain the observed large line broadening. 

1790 Appl. Phys. Lett, Vol. S1. No. 22,30 November 1987 


The observed line broadening of the spontaneous emis¬ 
sion is attributed to the shock-wave-induced exciton colli¬ 
sions. When a shock front propagated through the sample, 
all molecules behind it would gain a particle velocity u p via 
the intermolecular bonding, but excitons would not, since 
there was no such tight bonding between excitons or between 
exciton and molecules. Therefore, the emitting excitons in 
the shocked region all suffer from additional collisions with 
the entire array with a directional particle velocity u p . The 
collision frequency for each exciton is f c — a : u p N mo[ , where 
a is the Bohr radius of the exciton in GaSe which is 15 — 32 A, 
N mol is the molecular density in GaSe which is ~2x 10” 
cm -3 . The particle velocity at 17 kbar is about 2X10 4 
cms' 1 , Therefore, the collision broadening at 17 kbar is 
Av = 2f c /ir~ 2.5x 10 13 Hz or A/l — 33 nm which agrees 
well with the observed value. It should be noted that colli- 
sional broadening is homogeneous and should have a Lor¬ 
entzian line shape. 14 This is also clearly shown in Fig. 1. 
Since the collision frequency is proportional to u p which is 
related to the shock pressure P by the jump condition": 
P — P 0 = poU,u p , where U, is the shock velocity, p 0 and P 0 
are density and pressure in unshocked region. The shock¬ 
wave-induced collision broadening is expected to be propor¬ 
tional to the shock pressure. This is confirmed by our experi¬ 
mental data as shown in Fig. 3. The observed line broadening 
from different pairs of unshocked-shocked data plotted 
against the applied shock pressure displays a linear depen¬ 
dence. A pressure coefficient of line broadening for the spon¬ 
taneous emission of GaSe — 1.3 nm/kbar is obtained. 

The observation that the shocked spontaneous emission 
has the same integrated intensity as the unshocked one sug¬ 
gests that the elastic collision (phase-perturbing collision) 
plays a major role in our line broadening mechanism because 
inelastic collision (quenching collision) will decrease the 
number of emitting excitons and reduce the integrated inten¬ 
sity significantly. Note that the usual pressure broadening 
which is crucial for gaseous systems essentially comes from 
an increase in particle density and a decrease in mean free 
path due to high pressure. This latter effect cannot be signifi¬ 
cant in condensed matter because of its extremely small 
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FIG. 3 Observed wavelength broadening of the spontaneous emission spec¬ 
trum from GaSe vs shock pressure. 
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compressibility (typically p/p„~ 1.2 at /* = 200 kbar). The 
shocked-wave-induced exciton collision broadening is a new 
mechanism different from the usual pressure broadening 
and directly related to the unique nature of the shock waves 
and excitons in condensed matter. 
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